Zeitschrift fir 
angewandte Mathematik 
und Physik 


ZAM P 


Journal of Applied Mathematics and Physics 


Journal de Mathématiques et de Physique appliquées 


Editores: J.Ackeret EB. Baumann P.Scherrer E.S8tiefel F.Stissi W.Traupel H. Ziegler 
Redactor: R. Sanger 


INHALT —- CONTENTS —- SOMMAIRE 


Sears, W. R.: Rotating Stall in Axial Compressors. ... . et em ental chs Faery ts” eae 
SAuTScHI, W.: Uber den Fehler des Runge-Kutta-Verfahrens fiir die numerische Integration ge- 
wohnlicher Differentialgleichungen n-ter Ordnung . . ... +--+ + see eee ee eee 456 
PONDER, E.: Eine genaherte Behandlung des schweren symmetrischen Kreisels in nicht-Eulerschen 
moordinaten. 69. 965 2 ee eR Mears: eiien sah ioieiim lcs & \shisy foes 462 


Kurze Mitteilungen - Brief Reports - Communications bréves 


AntEeNneEN, K.: Uber den Mechanismus einiger Elektrete . . 6 ese ee ee ee et te es 478 
PrurpacuEr, H. R., und SANcEeR, R.: Mechanismus der Vereisung unterkihlter Wassertropfen 
durch disperse Keimsubstanzen iit tee | epee Pe EU Ress Teed mh e: Jet eT psc aicawel Pls 485 


Rora-Dresmnvtes, E.: Zur Berechnung der Geschossabweichung unter dem Einfluss eines Seiten- 


LCS eR ee ois) ah Siye pee e ll om ee Oe) Ns Operate A ats Poe mec eae! MEE Ree ae 494 
RUTISHAUSER, HL: prema tgs zur numerischen ae veweimbched Differentialgleichungen 

| (Gyan: "a fle Good a0 ght aie te gue aC ai eC Nd oa 497 
Sarro, H.: Torsion of a ereilan Shaft Pressfitted with a Disc’. . 9. . «© - ss ev a iets Mate? Ye 498 

Varia — Miscellaneous — Divers 

Herbsttagung der Schweizerischen Le hamranegs Gesellschaft vom 24. bis 25. September 1955 in 

EMITING weeipcpe sien en sl phe ee PRE Mall IAs Ite) Te 1k al saa ira it te pil wel te Niue de at) ¢ 503 
Buchbesprechungen — Book Reviews — Notices bibliographiques .. +++ +++ +2 ese 506 


Pag. 429—508 


Fasc. 6 25. 11. 55 | 


f\( Birkhauser Verlag Basel, Schweiz 


ee “4 a 5 


Zeitschrift fir angewandte Mathematik und Physik ~ 
Journal of Applied Mathematics and Physics 
Journal de Mathématiques et de Physique appliquées 
ZAMP 


Redactor: Prof. Dr. R. Sanger, Postfach Ziirich 23 (Schweiz) ‘ ‘ 


ZAMP erscheint regelmassig alle zwei Monate. Redaktionsschluss: spatestens 8 Wochen vor Erscheinungs 
_ termin. Manuskripte und alle die Redaktion betreffenden Zuschriften sind an den Redaktor zu richten. 
i Manuskripte fiir Kurze Mitteilungen k6nnen, falls die Autoren auf Zustellung eines Korrekturabzuges ver- 

__ zichten, 6 Wochen vor Erscheinungstermin eingereicht werden. Die Beitrage werden in deutscher, eng- ; 
He lischer, franzdsischer oder italienischer Sprache publiziert, doch soll jeder Originalarbeit und Kurzen 
t. _ Mitieilung eine kurze Zusammenfassung in einer andern als der Originalsprache beigegeben werden. 
et _ Abbildungsvorlagen sollen reproduktionsfertig mit ReduktionsmaBstab versehen eingeliefert werden, ; 
gegen soll die Beschriftung nur mit Bleistift, am besten auf einem lose vorgeklebten, durchsichtigen 
ve Papier, ausgeftihrt werden. Autoren von Originalarbeiten mit mehr als 16 Seiten Umfang haben von de 
17. Seite an einen Beitrag an die Satzkosten von Fr. 20.- pro Seite und von der 25. Seite an einen solchen 
£ 


von Fr. 35.— pro Seite zu entrichten. Die Verfasser von Zusammenfassenden Berichten, Originalarbeiten 
Kurzen Mitteilungen erhalten 50 Gratisseparata ohne Umschlag. 


ZAMP is published regularly every two months. All papers must be in the hands of the editor at least” 
1; 8 weeks before publication of the issue. Manuscripts and all other correspondence concerning contributio 
if to the periodical should be addressed to the Editor. Manuscripts of Brief Reports can be submitted 6 weeks 
2 fore publication of the issue, if the authors dispense with proof-reading. The contributions will be pub 
in English, German, French, or Italian; to each Original Paper and Brief Report should be added a bri 
summary in a language other than the original one. Drawings and photographs for illustrations shoul 
be sent ready for reproduction with indication of the proposed reduction; lettering on the illustrations should - 
be entered in pencil, preferably on a loosely attached transparent sheet. Authors of Original Papers exceed- 
ing 16 pages are requested to pay a contribution of Fr. 20.— per page for pages 17-24, and of Fr. 35.— per 
page for pages 25 and over. The authors of Survey Articles, Original Papers, and Brief Reports receive 50° 
free separata prints without cover. * 


ZAMP parait réguliérement tous les deux mois. Derniers délais pour présenter les manuscrits: deux mois 
avant la parution de chaque numéro. Tous les manuscrits et toute correspondance relatifs a la rédaction sont 
a adresser au rédacteur. Les manuscrits pour des communications bréves peuvent 6tre présentés six semaines 
avant la parution, si les auteurs renoncent A contréler les épreuves. Les articles seront publiés en langues 
francaise, anglaise, allemande ou italienne; cependant, chaque article original et toute communication bréve 
devront tre accompagnés d’un résumé en une autre langue. Les illustrations doivent étre prétes 4 la repro- 
duction et il faudra y indiquer l’échelle de réduction. La légende devra étre inscrite au crayon, de préférence 
sur papier transparent, légérement collé A l’illustration. Pour articles originaux qui dépassent 16 pages, 
une contribution de Fr. 20.— par page doit étre payée pour les pages 17-24 ; A partir de page 25 la contribution 
est de Fr. 35.- par page. Les auteurs de comptes rendus, ainsi que ceux d’articles originaux et de communi- 
cations bréves recevront 50 tirés A part sans couverture de leur article % 


en eee Sa ek ee : 
Jahresabonnement (6 Nummern): Inland Fr. 38.-, Ausland sFr. 48.—, Einzelnummer Fr. 10.— 


Annual subscription (6 issues): Switzerland Fr. 38.-, foreign Fr. 48.-, single copy Fr. 10.— 
Abonnement annuel (6 numéros): Suisse Fr. 88.-, étranger Fr. 48.~, le numéro Fr. 10.- 


Postabonnement — Post-Office Subscription —- Abonnement par office postal 
Belgien, Luxemburg Danemark Deutschland Finnland Holland 7 
Italien Matrokko Norwegen Portugal Schweden : 
Inserate/Advertisements/Annonces: 1/, Seite/page: Fr. 165.—, 1/) Seite/page: Fr. 88.—, 1/, Seite/page: Fr. 50.- 
a 
Alle Rechte, einschliesslich der Ubersetzung und Reproduktion auf photostatischem Wege oder durch 


Mikrofilm, vorbehalten. - All rights, incl. translations and photostatic and microfilm reproduction, re= 
served. - Tous droits, y compris traductions et reproductions photostatiques et par microfilm, réservés. 


Copyright 1955 by Birkhauser Verlag, Basel 10, Switzerland 


Printed in Switzerland 


Rotating Stall in Axial Compressors’) 


By Wri11am R. Sears, Ithaca, New York?) 


Introduction 


. . Experiments carried out in axial-flow compressors by several investigators 
mt]; (2},13]*) have disclosed that violently asymmetrical flow patterns occur. 


at) a ta eN 


in such machines when their flow rate is so restricted as to cause blade stalling. — Q 
“Typically, certain portions of the annulus of blades appear to be stalled while ©) 
“others remain unstalled; the stalled portions exhibit reduced airflow and the er 


_roughness of flow characteristic of stalling. But what seems most important is 
the fact that these patterns of stalled and unstalled flow do not remain fixed — 
to either the stator or the rotor blades, but rotate steadily in-the direction of © 
_rotor rotation at a lower speed than that of the rotor blades. The result is that. 
both rotor and stator blades are subjected to violent periodic aerodynamic 
loads, since they find themselves alternatively installed and unstalled flow... 


Be: “The experiments mentioned above [1], [2], [3] were carried out with the aid SS " : 
of hot-wire anemometers. One wire attached to a stator row is sufficient to — 


disclose that a periodic flow occurs; several wires, spaced circumferentially in 
ba stator row provide, in addition, information concerning the rate of rotation 
_ of the asymmetrical flow pattern. Patterns having one stalled region have been 
_ observed, and also patterns having several, equally spaced regions. The regions 
of stalled flow may extend across the annulus from blade root to tip (Figure 1), 
or they may lie near the blade roots (Figure 2) or near the blade tips. 
_. There are, in fact, interesting differences between the results obtained by 
_ the various experimenters. What is common to them is principally the fact 
F already stated above: that stalling of axial-flow compressors is accompanied 
_ by the appearance of steadily-rotating periodic stall patterns. This conclusion 
has been so generally reached that it may even be surmised that all axial-flow . 
~ compressors and pumps behave in this way. If this is correct, it may provide 
an explanation of numerous fatigue failures of blades, which have previously 
“been blamed on flutter or stall flutter, or on other causes. 
It should also be made clear that this phenomenon is not the familiar one 
_ known as ‘surge’, which involves periodic variation of the rate of flow through 


? 


1). This research was partially supported by the United States Air-Force under Contract 
No. AF 33 (038)-21406, monitored by the Office of Scientific Research, Air Research and Develop- - 
ment Command. — ; 
2) Graduate School of Aeronautical Engineering, Cornell University. 
- 3) Numbers in brackets refer to References, page 454. 


Figure 1 


Depiction of rotating-stall pattern [2]. Three-stage compressor. 


em Figure 2 
. bays Depiction of rotating-stall pattern [2]. Three-stage compressor. Slightly higher flow than for Figure 1. 
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the machine. This distinction has already been pointed out clearly [1], [2], but 


“still seems to require emphasis. The phenomenon of rotating stall is basically 
one involving steady flow through the machine; only its distribution around the 


annulus is nonuniform. The flow, in fact, becomes a truly steady flow when 


viewed from a coordinate system rotating with the stall pattern—a fact that 


will be exploited farther below. To be sure, the two phenomena, rotating stall 


and surge, may occur together, and it appears that the non-uniqueness of the 
-rotating-stall type of flow for any certain flow rate may be one of the chief 
_ sources of the instability that we call surge [2]. 


The present paper is a review of several attempts to explain the phenomenon 
of rotating stall on a theoretical basis. As will be seen, the attempts have not 
been entirely successful to date, especially because they concern a type of 
flow — unsteady, stalled flow — about which there is insufficient knowledge. On 
the other hand, they are not completely unsuccessful, since they show, at least 
qualitatively, that flow patterns of the rotating-stall type may occur in axi- 


symmetrical rigid blade rows with steady, uniform inflow to the machine. This 


conclusion is of some importance. It seems clear that stall flutter, for example, 
could produce periodic flow patterns that would rotate around the blade 


annulus if the flutter frequency were not the same as, or a multiple of, the 


rotational frequency. But the occurrence of stall flutter must depend intimately 
on the elastic properties of the blades. Thus it is important to know whether a 


_ similar phenomenon can occur purely aerodynamically, in a row of rigid blades. 


mS 


Notations 
constant in equation (35) ; 
V/U (Figure 3) ; 
constant in equation (35); 
oc; lift coefficient of blade; c,{«] denotes the lift coefficient as function 

of «, while c,(y) represents the same quantity as function of 4; i.e, 
cyla(y)] = ex(y); 

C constant in equation (35) ; 

C,,C, constants in equations (8) to (12 

h hy) = AH [o U*; 

h, Fourier coefficient of A(¥) ; 

H H(x, y) = total pressure ; 

AH AH(y) = increase of H across blade row at y; 

k circumferential speed of blades relative to flow pattern, divided by 


axial flow speed U; positive in negative y direction ; 
K,, K, constants in equations (38) and (42) ; 
L circumference of blade wheel; 
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(local) slope of 4 against 6; in channel relations; 

coordinate normal to undisturbed stream direction [Figure 3 and 
equation (3)]; 

summation index of Fourier series; 

p(x, y) = static pressure ; 

Ap(y) = static pressure rise across blade row; 

real part; 

time; 

u(x, y) = perturbation axial velocity component; 

Fourier coefficient of u(y) ; 

u(0, ¥); 

undisturbed axial velocity component ; 

v(x, y) = perturbation circumferential velocity component; 

v(— 0, 9); 

BOSH); 

bU undisturbed circumferential velocity component; 

axial coordinate (Figure 3); 

circumferential coordinate (Figure 3); 

a(v) = (B; + B)/2; mean angle of attack; 

Fourier coefficient of «(y) ; 

B.(y) inlet angle relative to blade (Figure 4); 

mean (unperturbed) value of f;; 

B.(y) outlet angle relative to blade (Figure 4); 

Fourier coefficient of 6,(y); 
angle of phase lag of h behind £;; 
angle of phase lag of c, behind «; 
variable of integration; 

flow coefficient [equation (36)]; 
density of fluid; 

solidity of blade row (= chord/gap) 
w(x, ¥) = stream function; 

Q(x, y) = vorticity. 


en SF 
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(local) slope of the curve c,[a] (per radian) ; 
. 
: 
4 


’ 


Subscripts, superscripts, etc. 


value at x = — oo; 
value at x = + co; 
value given by the uniform-turning solution; 
rotational part of (~); 
irrotational part of (~); 

variable part of ( ), defined by ( ) 


I 
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Flow Through a Blade Row 


All of the theoretical investigations of rotating stall, to date, are based on 
four simplifying assumptions; viz., 

(1) The compressor is replaced by a single blade row. 

(2) Assuming that the blade length, i.e., the width of the blade annulus, is 
small compared to the mean radius, the blade row is replaced by a two- 
dimensional cascade of blades. 

(3) The flow is supposed to consist of a small perturbation of a uniform stream 

' (4) Viscosity and compressibility of the fluid are neglected. 

In all but one study [5], the further approximation has been made that: 

(5) The blade row consists of infinitely many blades; i.e., the row is replaced 
by a discontinuity surface at which energy, due to rotation of the blades, 
is introduced to the flow. 


Figure 3 


Diagram showing coordinate system and notation. 


In the lone exception, which is an unpublished investigation [5], a row of 
discrete blades was first treated. To render the results practically useful, 
however, the same approximation of vanishing chords and blade spacing was 
later introduced, so that the work finally fell into the same category. 

The two-dimensional incompressible-flow problem resulting from these 
assumptions is sketched in Figure 3. The discontinuity surface lies along the 


ZAMP VI1/28 


434 WILLIAM R. SEARS ZAMP 


y-axis. It will now be assumed that the flow pattern is a steadily rotating one, 
as has been described above, and the coordinate system will be defined to be 
fixed relative to this moving pattern. The flow is thereby rendered steady, 
with a consequent simplification of the analysis. For example, BERNOULLI’S 
equation can now be used to describe the flow either upstream or downstream 
of the discontinuity sheet. It should be recognized that the analysis applies 
equally well to a stator or a rotor row, since there is no distinction between 
these from the point of view of the coordinate system employed. 

Other notation to be employed in the analysis is indicated in Figure 3 and 
in Notations, page 431. So far, the problem is formulated exactly as in Reference 
[4] except that the swirl angle relative to the moving coordinate system, tan! b, 
is not assumed to be small. The appropriate, linearized Bernoulli equation is, 
for 2.=..0;, 


pteUut+eVuv=f_,, (1) 


where u and v are the perturbation velocity components and #_ ,, is the pressure 
far upstream where conditions are uniform. 

Let the total-pressure rise through the blade row be AH(y). Simply to 
provide a dimensionless representation of this quantity, let AH(y)/(@U2) be de- 
noted by A(y). The Bernoulli equation for the downstream region is then the 
statement that p+ @Uu-+e@Vv is constant and equal to p_. + AH along 


any streamline. Within the scope of a small-perturbation theory, however, this — 


statement can be revised as if the streamlines lay in the direction of the 
undisturbed flow. Thus, for «> 0 , 


pP+oUureVv=p_~+oU*h(y—b x) =H(x, y). (2) 


7 


In Reference [4] the vorticity in the flow downstream of the actuator was 
calculated by consideration of the shedding of vortices by the blades as they 
pass through a variable flow pattern. Actually, this resort to airfoil theory is 
unnecessary, since for steady flow there is a simple relation between the vorti- 
city 2 and the gradient of the total pressure H. This relation is, for the present 
application [10] . 


, 


stipe G5: AGT Le) (3) 


On 


where 0H/dn denotes the derivative of H in the direction normal to the undis- 
turbed streamlines; viz.., 


0H OH U OH V 


on ~ dy VuRp Ve Ox orp ~ CUM (y—b)VUR+Y?. 
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Therefore the vorticity is 

Q(x, y) =—Uh' (y—b x). (4) 
The boundary conditions far upstream are the vanishing of all perturbation 


quantities. The boundary condition far downstream is the condition of uniform 
pressure and parallel flow. The condition of parallel flow is 


V+u V v u 
meena ni iz (1 ei 7) = constant (5) 
or 
Uv U 
7 > yu = Constant . - (6) 


It should be noticed that the direction of this parallel flow is not specified; it 
differs, in fact, from the direction of incoming flow. Nevertheless, this condition 
is Sufficient to establish conditions at x = oo, except for constant terms not 
required in the analysis, for equation (6) leads to 


dv Ou Ou Ou 
areas Ri nee bp (7) 
at % = 00. 
With Q given by equation (4), it is easily verified that the solution of this 
differential equation for u is 
Ugg = U (1+ b%)-*h(y— bx) + Cy, ea) 


and therefore 
v,, =U b(1+b?)1h(y— bx) + Cy, (9) 


_ These values, u,, and v,,, are the asymptotic values of and v far downstream, 
while C, and C, are constants involved in the uniform-turning part of the 
solution, which is discussed below. It is easily verified, by substitution in 
equation (2), that these yield a constant value for px. 


The Uniform-Turning Solution 


Having established the flow conditions far upstream and far downstream, it 
is now possible to identify a solution consisting entirely of constants, which 
represents uniform turning of the flow through the blade row. Let these values 
be denoted by the subscript c. They are 
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Fon 4< 0: | 
Uu-=0=4%,, Pe = Poo 
Borin >.0: 
1 b 
Ui Oe as Uh,+C,, ams py 3 Uh,+C,, (10) 
c 1 
Po Pao aon Ceres U? he = Poo 
or , 
0U? hy= ~,, — Px + oVu, 
and | 


=e 


The Boundary-Value Problem 


The uniform-turning solution is a trivial one so far_as the phenomenon of 
rotating stall is concerned. To proceed to more interesting solutions it is neces- 
sary to work out the relations that exist between the flow variables both 
upstream and downstream of the blade row. These are, of course, different in 
the two regions because the flow upstream is irrotational while that downstream 
is rotational. The two regions are related at the blade row (x = 0) by the con- 
tinuity of the axial flow component and by the characteristics of the blades 
themselves, which will be discussed later. Let the velocity components just 
upstream of the blade row be called u(y) and v_(y). Then the value of « just 
downstream of the row must also be u(y), while the tangential component will, 
in general, jump to a different value v aly}, 

In terms of the familiar plane stream function p(x, y) whose derivatives 
are y, = u and y,=—v, the boundary-value problems of the upstream and 
downstream regions are, for x <0: 


ras Wear NAGS ear oan | 


(11) 
p(—°o, y) = _., =constant, u(—0, v)=(y), v(—0O, yy =v ly), | 
Bora 0: 

Veyp=-2=UR' (y—b x), 
1 
u(oo, y) = EERE Uh(y—bx)+C,, 
b 
v(co, y) = ibe Uh(y—bx) + C. Ma) 


Pico, y) =~ _..- 0 UC,— oVC,, 
u(+0, y) = u(y), v(+0, y) = v4(y) . 
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For convenience, let all these quantities be considered as made up of the 
-constant-turning parts, equations (10), and variable parts denoted by the 
sign (_); i.e. p= y,+ y, etc. Then, since the former consist of undisturbed 
flow upstream, equations (11) are the equations for y, u, v, and , except for 
the boundary condition p(— ce, y) = 0. Equations (12), on the other hand, are 
sreplaced, for x > 0, by: 


u= it; uloo, 9) = zap UA (y—ba) +, 
il ~ 
u(oo, y) = iaetEs Uh(y— bx), eo 
~ b ~ 
0(00, ¥) = 44 ga Uh(y— 5), 


b(co, y) =0, u(+0, y) = u(y), 
v(+0, y) =v,(y) — % = 24(y)- 


a 


The solution of equations (11) can be written down immediately in the form 
of definite-integral relations between u(y) and v_(y). One can consider these 
as resulting from a distribution of sinks along the y-axis, of strength 2 Uo(¥) Or 
from a distribution of vortices of strength 2v_(y). The results, which are 
strictly equivalent, are 


ry d il ae n) a1 
Uo(¥) = 2 / rae and. v_(y)=—< if ah (14) 


where the integrals are Cauchy principal values. 

The solution of equations (13) can be constructed by the same technique 
after a particular solution of the Poisson equation is subtracted out. A suitable 
rotational flow for this purpose is the parallel shear flow far downstream, given 
by i1(co, ¥), 0(co, y) in equations (13). In other words, let 


° 1 _ 
“=u 44? with wv? = —, Uh(y— dx) (15) 


and 
p= oO) 4 yl?) swith,...v) = ‘a Uh(y— 6x). (16) 


~~ 
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Then, with y = y") + y®, one has for the downstream region, 


w(+0, 9) = toy) — 40, 9) = ol¥) — pe U ALD) (17) 


b ~ 


. 
V2y—0,  ul(o0, ») = 0 = v"(00, y) 
v0, y) =i — U0, 9) = 0. — Sa A) 


: 


The relations analogous to equations (14) are then 


eee aa Lf valah— Uta Eos 
u(y) — ight Ony) = — = / _ y= ¥ dy 
and (18) 
i b ~ 1 fugly) — UR(m (1 + b2) 
Bath aeae Ui Lf Soha UND gy, 


Thus, the velocity components upstream of the blade row are definitely 
related to one another, while those downstream are related in a way that 
involves the total-pressure-rise function h(y). 

All the small-perturbation analyses of the rotating-stall phenomenon are 
basically the same up to this point, although some of the authors have not 
made use of the steady-flow coordinate system, and have had to look for those 
special unsteady solutions that do not grow or decay as they propagate. At this 


point one has just two relations, equations (14) and (18), between four unknown 
functions, 


u(y), v(y), d,(y), and Ahly). 


The other two relations needed may be called the ‘blade-characteristic’ rela- 
tions; it is in the choice of these that our ignorance of the precise nature of 
stalled flow becomes apparent and in which the various analyses of the phenom- 
enon diverge from one another. 

Before proceeding to discuss the various assumptions that have been made 
regarding blade characteristics, it may be of interest to consider briefly the 
physical and mathematical natures of the solutions sought. The fundamental 
unknown function is the pressure rise hy), which describes the loading of the 
blades. The equations derived above state the relations that exist between this 
loading and the induced velocities Uo(¥), v_(y), and v,(y). The blade-charac- 
teristic relation will tell how the loading is affected by variations of the velocity 


and flow direction at the blade row; i.e., ultimately, how hy) is affected by the 


me 
e 
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_ same induced velocity components. Thus, the problem is strictly a homoge- 
neous one: What nonuniform blade loading function induces such velocity 
“components as just to support itself ? 
There is an analogous problem in the realm of lifting-line wing theory: What 
‘spanwise circulation distribution induces such a downwash distribution as just 
to support itself? It is easy to show (and seems physically obvious) that for 
unstalled wings PRANDTL’s equation admits no trivial solution to this problem. 
If the slope of the section lift curve is negative, however, the homogeneous 
_ Prandtl equation has infinitely many eigen-solutions, but only for eigen-values 
of the slope [11]. One is not surprised to find that the rotating-stall problem 
: has analogous solutions, which will be discussed farther below. 


Blade-Characteristic Relations of Airfoil Type 


In References [4] and [5], the two additional relations needed to connect the 
upstream and downstream regions of flow were carried over from airfoil theory. 
It was envisaged that, even in conditions of partial stall, the most important 
force acting on a blade would be the lift, or more precisely, the component 
arising from circulation. This is the force that, in a stationary cascade, would 
produce a pure deflection of the flow without gain or loss of total pressure; in 
a moving cascade, of course, it results in a change of total pressure. Taking 
into account this force alone, the first blade-characteristic relation was written 
by equating the rise of total head to the axial force exerted by the moving 
blade. Since the circulation per unit length of the discontinuity sheet is v_ — v,, 
this relation is 

AH(y) = @ RU [v_(y) — v,(y)] | 
or | (19) 


hy) = fea) = eso) 


where k U denotes the velocity of the blades relative to the moving coordinate 
system. : 

The second relation, according to this formulation of the problem, is that 
connecting the blade circulation v_ — v, to the blade angle of attack «: 


-(y) —0(9) = 4 o UV 1+ 0+ A? ela), (20) 


where o denotes the solidity of the blade row and c,[«] is the lift-coefficient 
curve of the blades, which is nonlinear, in general. According to the theory of 
airfoils in cascades, and neglecting unsteady-flow effects on the ground that the 
wave lengths of flow patterns are large compared to the blade chord (already 


: ' 
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neglected), one should use here the mean angle of attack ; 
i : 
a= Abe PBs) ; 


It is easy to compute this angle as the sum of arc tangents and then to linearize 


2 aurwer, 


\ 


F Ut+u,” 
Ga 
kUebU ev. 4 a 
Utuy 


Figure 4 
Diagram showing inlet and outlet angles, B; and £,. 


the result in accordance with the small-perturbation hypothesis. The result 
gives « in the form «,-+ « where «, corresponds to the uniform turning solution 
discussed above. The variable part, «, is found to be 


1/2 (v_+ 04) — (k+ b) uy 


Axe : ay 
Uli + + &)*) pe 
The second relation may now be written in terms of perturbation quantities: 
= 1 ae = ieee 
vy) — H,(y) => oUV1+ (b+ 8)? G fa], (22) 


where ¢,|«| is, of course, the lift-coefficient function measured from the constant- 
turning point cy, = ¢[a,] (Figure 5). 


Linearized lift-coefficient curve. Phase lag 


It seems most consistent with the small-perturbation assumption to linear- 


ize the curve at this point, i.e., to replace the curve ¢,[a] by its tangent, and 


to write 
ela(v)] =ma(y) (m= constant) . (23) 


This is the type of approximation made in Reference [4] and included in a 
more general framework in Reference [7]. In both cases, moreover, the phenom- 
enon of time lag was also included in an approximate way. It is well known 
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Figure 5 


Diagram defining lift-coefficient functions ca] and [a]. 


that the circulation about an airfoil lags behind its incidence by an interval 


which is greater at large angles of attack than at small. This time lag is due 


primarily to the time required for the viscous processes of separation and 
vortex shedding after a change of incidence, and may be considerably larger 


than the subsequent time lag in the build-up of circulation according to the 
theory of airfoils in unsteady motion. It has been studied by MENDELSON [6], 


who concluded that, approximately, the phase lag for any airfoil in oscillatory 


motion depends on the mean incidence and not on the frequency of oscillation. 

This interesting result presumably cannot apply to all frequencies, but was 
proposed as an approximation for the range of frequencies involved in many 
cases of stall flutter. In the present application it states that, if the angle-of- 


attack perturbation is sinusoidal: 


a(V) = ER {on iad i si : ey 


where «,, ”, and / are constants, / being the circumference of the blade wheel, 


then the corresponding lift coefficient is 
¢,(y) = ayn Ro, EE RN ob (25) 


where the phase-lag angle A depends upon «, but not upon «, or 7. (Here the 
symbol “Y has been used to denote the real part of a complex quantity. Inwhat 
follows, this symbol will be omitted, and real parts always implied.) 

To utilize these relations, suppose /(y) to be expanded in a Fourier series. 
It must have a period equal to /. Thus, 


oO 


h(y) ae hn e2nin yt (26) 


n=1 


ne 
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The corresponding induced angle of attack «(y) then has the form 
ys 
; 


= Yay ein, (27) 
’ n=1 

and equation (25) states that | 
co 

Ey) =m >: a, erm (28), 

n=1 A 


But the Fourier coefficients «, are easily found by means of the equations 
derived above, viz., equation (21) for « together with equations (14) and (18). 
In particular, one finds that 


b ‘ i eee ne 
OD) a Dey ha rere; PG) Ne ee 5) J gy 
a) (29) 
b+1 = 2niny/l ; 
a5 hye 
and 
_ 1 fen — oly) Ue URN PB ae 
U 1 -i(n) dn 2rin y/l (30) 
=5 =f a +O DT he e 
— |i oe 18) Qt fy 
a cal Te | | She e 
so that 
gs 2) hal fae 


a(y) = ie Oth? 1 1+b2 ils Odes aloes ezting/L (31) 


The coefficient of e**'"*/’ in this expression is the Fourier coefficient «, 
of equation (27); thus ¢,(y) follows from equation (28), and equation (22) 


becomes 
co 


D> hn erring lt _ ; mak Vi-+ (b 7 Ry? FS, et Ran(y/l) +4] (32) 


n=1 


This is the complete equation of the problem, according to the present 
approximations. It is the statement that the loading h(y ) is supported by its 
own induced velocities. It is, as was mentioned above, a homogeneous equation, 
and can be satisfied in a nontrivial way only when the successive Fourier 
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‘constants on the right- and left-hand sides are equal; 1.e., 

1 Seecee ORE eT 3 

hy = mo kV1+ (+h)? a, e'4 (33) 


with «, taken from equation (31). This requires both 


k2b uae 
appr t2k+ b= 775s tana 


1+ 


and (34) 


2 ——— oe 
“ = oe V1 | (b+ k)? cosA. 


The simultaneous solutions of these two equations, for various assumed 
values of A, are plotted in Figure 6. The most convenient independent variable 


Figure 6 
Rotating-stall solutions based on blade-characteristic relations of airfoil type. The abscissa, k + D, 


is the tangent of the mean inlet angle Ba Where double-valued functions occur, corresponding 
values of & and —mo/4 are to be read from branches showing arrowheads in the same direction. 
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: 


is k + b, which is the tangent of the mean (undisturbed) angle of inflow relative : 
to the moving blade; i.e., the mean value of #; in Figure 4, say f;. This is ae 
parameter of clear geometrical meaning, easily estimated from design data for 
any given blade row in a turbo-machine at any known flow rate. Equations (34) 
state that, for any given value of this parameter and for any known value of A, 
rotating stall can occur only for two special values of m o—both negative—and 
that the stall pattern must rotate at certain corresponding speeds & U relative 
to the blades. 

The solutions plotted in Figure 6 are a generalization of those presented in 
' Reference [4], where a restrictive assumption of small swirl (6 << 1) was made. 
It is important to notice that no information is obtained from the solution 
concerning the amplitude, shape, or wave length of the periodic solution h(y). 
This is a result of the linearization of the problem, including the blade-character- 
istic curve. According to this result, any periodic function h(y) ‘supports itself 
by means of its own induced angles, provided only that k + 6, mo, k, and A 
satisfy equations (34). 

\ 


Nonlinear Lift-Coefficient Curve [5] 


In Reference [5], MARBLE attempted to avoid this situation by introducing 
a certain nonlinear lift-coefficient function in place of equations (23) and (25). 
This function, suggested by H. S. Tsren, is 


ala]=Aa—B(a—C - tre (35) 
where A, B, and C are constants and da/dt denotes the time rate of change of «. 
He then assumed h(y) to be a simple sine, and tried to find solutions of the 
rotating-stall problem in a manner analogous to what has been carried out 
above. Preliminary results were presented in Reference [5]. 

In general, the nonlinearity produces higher harmonics in the induced angle 
&, even when /i(y) is sinusoidal; thus the problem of satisfying the homogeneous 
equations is a much more complicated one. In other words, a mathematical 
difficulty arises in the determination of the Fourier coefficients of h(y), since 
they are not independent of one another. In this situation Professor MARBLE 
has concluded that it is more fruitful to use a different description of the 
nonlinear blade characteristic at stalling, which will be discussed below. 


Blade-Characteristic Relations of Channel Type 
In some other investigations of the same phenomenon, the concept of the 


blade row as a lattice of airfoils has been discarded. Instead, the row has been 
treated as an ensemble of channels. 


Se 
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3 Variable-Area Channels ({1] and [9}) 


- 


__ In Reference [1], Emmons suggested that the effect of partial stalling of the 
[ channel bétween two blades could be represented as a partial constriction of 
_the passage area through the channel. He defined a flow oat * for the 


passage such that 
U + u%) = const - x. (36) 


This was based on the assumption that the static pressure behind the blade 
row would remain uniform in spite of the asymmetric flow pattern, so that the 

only variable affecting the flow through the row at any point would be the 
degreee of constriction, i.e., the degree of boundary-layer separation. EMMons 
further assumed that the flow coefficient is determined by the inlet angle 6,(y) 
For the purposes of his investigation he linearized the expression for this angle, 


viz., bak os S 
p:= 1+ (6+ R)? Fes A tot) 


and assumed a linear relationship between x and f; in steady flow. For unsteady 
flow, he assumed that x increases at a rate proportional to its deviation from 

this steady-flow value, for any given es this is equivalent to assuming a 
certain type of lag between B, and %. The result of these assumptions is to 
provide a blade-characteristic function in the form 


u'o(y) = Ky u% — Ky v=, (38) 


where K, and K, are real constants involving the various constants of propor- 
tionality of Emmons’ theory as well as # and 0. 

According to this theory, it will be noted, the flow pattern upstream of the 
blade row is independent of the flow downstream. Equations (18) are not 
needed; the only unknowns are u(y) and v_(y), and they are determined by 
equations (14) and (38). Assume that m(y) is given as a Fourier series (real 
parts being implied, as before) ; then 


= dD, Mn ee AL (39) 
n=1 

eee d oo 
v_(y) << 1 ( sol EA = iS, or*™ y/i (40) 

2. oan, t n=1 

274 = Qniny/l 
U y) oy >” Le 2 ae (41) 
n=1 


so that equation (38) requires 
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Emmons has already noted [1] that this means that K, =0 for a flow 
pattern propagating without growing or subsiding. This amounts to a condition 
on the slope of the curve of x versus ,; it seems to be analogous to the condition 
on mo obtained in the earlier analysis [equation (34)]. He did not mention. 
that it also requires a special relationship between m and Kj; i.e., his theory 
states that the rotating-stall pattern must be purely sinusoidal with a wave 
length determined by the constants of his blade-characteristic relations and 
parameters k and 6. The nature of this relationship is such that the relative 
speed k U decreases with increasing n, i.e., with decreasing wave length of the 
pattern, all other parameters being the same. : 

This theory has been extended by STENNING [9]. STENNING’S analysis is 
basically the same as Emmons’, but an attempt is made, in introducing time 
lag, to distinguish between ‘boundary-layer lag’ and the time lag due to inertia’ 
of the fluid contained within any blade passage ; it is concluded that the latter 
is more important. The same assumption of uniform downstream static pressure > 
is made and is defended as an approximation for stall patterns whose wave 
length is not large compared to the circumferential blade spacing. 

STENNING’S results are similar to Emmons’. Again two relations are found 
for steadily rotating stall patterns, one of which is a condition on the slope of - 
the curve of flow coefficient x. The other (which is fully discussed by STENNING) — 
restricts the steady pattern to a pure sine curve and relates k, b, und l/n in such 
a way that sinusoidal patterns of long wave length rotate faster relative to the 
blades than those of shorter wave length, all other effects being the same. 
STENNING presents some experimental data, obtained by Professor Emmons, 
which seem to support this conclusion regarding the effect of wave length on 
speed of propagation. There are other experimental data, however (e.g., in 
Reference [2]), which do not show the same tendency. 


Channel Relations Used in Reference [7] 


Besides the ‘airfoil’ analysis already discussed above, Reference [7] gives 
a solution based on ‘channel’ equations different from Emmons’ and STEN- 
NING’S. Following a suggestion of Professor RANNIE, the author undertook to 
set up blade-characteristic relations that would resemble the results of steady- 
flow cascade tests. Such results are sketched in Figure 7, where it is indicated 
that (1) stalling of a stationary cascade is indicated by a sudden increase of 
total-head loss, — AH, plotted against inlet angle B,, and (2) the outlet angle p, 
is substantially constant, for high-solidity cascades, even in the stall. The first 
of these conclusions must be modified to apply to a moving cascade, of course. 


For this case the corresponding curve of AH versus 6; would have the character 
of Figure 7 (c). 
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These curves suggest the use of the following blade-characteristic relations: 


~ 


E- 7 By =<cOnstant, OF B. =0 (43) 


and, linearizing the curve of Figure 7(c) at the appropriate uniform-turning 


Figure 7 
Sketches showing (a) outlet angle and (b) total-pressure rise of typical stationary cascades and 
(c) corresponding total-pressure rise for rotating cascades. B; is the inlet angle. 


values, 
h=—M£8, (M=constant) . (44) 


The sign convention adopted here is such that M is positive for the stalled 
blade row. 
As a generalization of these relations, the type of time lag (phase lag) obser- 
_ ved by MENDELSON for airfoils was introduced at this point in Reference [7]. 
This was, of course, rather arbitrary. It is justified only by analogy with the 
airfoil result. The process of viscous separation and energy loss described by 
the channel relations is similar to the process of airfoil stalling, which MENDEL- 
son found to be approximated by a phase lag of this type. Furthermore, it 
seems to be the only type of lag that can be introduced without automatically 
restricting the flow pattern to a simple sine. Rather, the enormous generality 
of Fourier series is maintained; Le., solutions will be found for which 


ie (45) 
n=1 
and 
h(y) ig. Ms Gree ae Threl (46) 
n=1 


where 6 is a real constant. 
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Again the process of solution begins by finding the Fourier coefficients B,, 
corresponding to the /,, of equation (26). The linearized expression for f, is the 
same as equation (37) but with #, written in place of B; and @, in place of v_; 
thus the first channel relation, equation (43), states that 0, = (6+ k) uw, and 
6;, from equation (37), becomes 


= 1 v_—v 
Bly) = 1+ (b+ k)2 ae eh Gs +. (47) 


| 
| 


Now, an expression for v_— 9, involving up, is obtained from equations (14) 
and (18): 


? 
4 


~ 


0-9) —H(9) = — = [Holden Fi 4 t. Uf Alen ays 


a y—n 1406? 1+ 6? y—-n - 


ee —oo 


Another expression for the same quantity follows from equation (14) and (43): 


co 


- 1 Uo() ay, : | 
vty) — Bely)—= — Ff MMSD — 4b) rly) (49) 
One can obviously eliminate v_ — 3, from these two equations, obtaining 
1 fulndan BOM eae Lf h(n) an | | 
zp yang ~ OFF) Ho) = rae {— daly) +2 f AMA2), oy 
If #,, denote the Fourier coefficients of Uo(¥), So that 
Uo(V) aes tke ezrin ft 
equation (50) relates w,, to h,,: 
U a+b 
Uy = 1 + 62 2 :o o k h,, . (51) 


Now, making use of equation (49) again, one has a simple relation between u,, 
and the Fourier coefficients of v. — 3,, which in turn are related to fp, i 
equation (47). Finally, from equation (46) the homogeneous equation of the 
problem becomes 


fo. 3) 


2Zainy/i __ 1 1 5 ~ [2a HT 
J hae Me = Ot A] Sr ugeilene ven 


nm=1 


(52) 


, ; co 
ros M : t= (64+) a+b Sig Ce MO EO 
=1 


i} 
14+ (b+h)2 i+btk 1482 . 


i 
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_ This requires both 


_ ie 2 

4 ieot, 1+ (b+) (53) 
and 

”, ; ins 1 acho: he (54) 


b(b+ kh)? 4+ 2k+ 0° ; 


When equations (53) and (54) are both satisfied, any periodic flow pattern is in 
equilibrium and rotates steadily. The compatible values of M,k, and b+k 


i 
‘ 


Figure 8 
Rotating-stall solutions based on blade-characteristic relations of channel type. The abscissa, 
k + b,is the tangent of the mean inlet angle B,. 


are plotted in Figure 8 for several values of the phase-lag angle 6. Again the 
independent variable is 6 + k, the tangent of the mean (undisturbed) angle of 
inflow relative to the blades, tan f;. All the interesting values of M are positive, 


_ ZAMP VI/29 
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which means that the blade row must be so fully stalled as to have a negative 
characteristic slope, as defined in Figure 7(c). Then rotating stall can occur at 
a certain relative speed k, with or without lag. ; 
= 
Marble’s Discontinuous Blade Characteristic . 
In Reference [8], Professor MARBLE has employed a nonlinear blade char- 
acteristic of a most interesting type, namely, a relation exhibiting a disconti- 
nuity at the stall. That such characteristics should be used was suggested in 
Reference [4] but was not carried out. MARBLE’s first blade-characteristic 
relation is a generalization of equation (43); viz., . 
By = const - B;. (55) 
His second relation is expressed in terms of the static-pressure rise across 
the blade row, , 
Ap = p(+ 0,y) — p(— 0, y) - 


It is his assumption that this pressure rise varies approximately linearly with 
the inlet angle f; below the stall but drops to zero when stalling occurs. This 
relation is sketched in Figure 9. 


Ap 


Ss 


Figure 9 
Sketch showing MArsie’s discontinuous blade-characteristic function [8]. Ap denotes static- 
pressure rise and B; the inlet angle. 


For convenience in applying this condition, MARBLE works with the depend- 
ent variable (x, y) instead of the velocity components used in the previous 
investigations. Since the small-perturbation approximation is employed, the 
pressure perturbation satisfies LAPLACE’s equation both upstream and down- 
stream of the blade row. MARBLE shows that the angles £, and Bo can also be 


tS 


a 
See 
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"expressed quite simply in terms of the pressure perturbation and its harmonic 
_ conjugate. Thus, the boundary-value problem of the two-dimensional blade row 
is reduced to finding a harmonic function whose real and imaginary parts, 
along the axis x = 0, satisfy certain relations derived from equation (55) and 
the characteristic sketched in Figure 9. 
_ The results of this analysis will not be reproduced here. They consist of the 
propagation speed ratio k and the percentage of circumference stalled, both as 
_ functions of k + 0b, the tangent of the mean inlet angle. MARBLE finds that & is 
_ insensitive to variations of k + b for values of k + 6 near 1, which is a typical 
value for axial-flow compressors. It is interesting to notice that in one parti- 
cular case that he studies in detail, in which the constant of equation (55) is 
taken equal to zero, so that this relation coincides with equation (43) of the 
earlier, SEARS analysis, his curve of & versus k + 0 is identical with the one 
_ plotted in Figure 8 for 6 = 0. Presumably the effects of lag could be introduced 
into MaRBLE’s analysis, but this has not been done. 
The most interesting result in MARBLE’s work is the percentage of circum- 
ference stalled, or the ratio of stalled to unstalled length in the periodic flow 
pattern. This increases rapidly from 0 to 1-0 as k + d is increased above the 
_ value where partial stalling first appears. The slope of this increase depends on 
the magnitude of the pressure rise 4p across the unstalled portions but seems 
to be relatively insensitive to the slope of the unstalled portion of the blade 
_ characteristic. 


Comparison with Experiment 


Practically the only comparison that can be made with experimental results 
at present is the comparison of measured and predicted rates of stall propaga- 
tion. Such a comparison is attempted in Figure 10, where the curves of k versus 
k + b are collected from Figures 6 and 8. The experimental points in this 
figure are taken from the report of IurA and RANNIE [2]. They represent the 
conditions of the rotor and stator rows of their experimental compressor at 
the onset of rotating stall. It is not known whether the rotor or stator rows 
were responsible for the stalling. To be more precise, it must be admitted that 
the present single-row theories are only roughly applicable to IuRA’s and 
Rannie’s three-stage machine, and the attempted comparison is based on the 
assumptions (1) that either stator or rotor rows (but not both) must be respon- 
sible for the onset of rotating stall and (2) that the presence of other rows does 
not seriously affect the results, as compared to those of the idealized single- 
row model. 

The dashed lines emanating from the experimental points in Figure 10 
show the approximate variation of k with k + b, for rotors and stators, respec- 
tively, for flow quantities smaller than that at which rotating stall commences. 


7" 
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Figure 10 
Comparison of measured rates of stall rotation [2] with theoretical values based on airfoil (above) 
and channel (below) types of blade characteristics, from Figures 6 and 8. The abscissa, k + 3, is 
the tangent of the mean inlet angle B,. The experimental points and dashed lines are from Reference 
[2], computed for rotors and stators, respectively. 


It is seen that the experimental points fall within the range of the theoretical 
diagrams based on either ‘airfoil’ or ‘channel’-type blade characteristics, 
However, the point denoting rotor conditions requires either negative lag or 
excessive positive lag if the ‘channel’ theory is adopted, while the point 
denoting stator conditions requires a rather large value of the lag (compared 


+ chee me Oe 


a 


, a ; ; 
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is to MENDELSON’s observations) in the ‘airfoil’ theory. There is an interesting 
_ coincidence between the experimental curve for rotor conditions (the upper | 
_ dashed line) and the theoretical curve for constant phase lag A in Figure 10a. 
_ Nevertheless, it is clear that the available comparison with experimental results 


_ does not lead to any general conclusions. 
my 


Discussion and Conclusions 


oc 


In this paper, several attempts to explain theoretically the phenomenon of 
rotating stall have been described. These are all based on the assumptions of a 
single row of closely spaced blades and the small perturbation type of flow. 
They are, therefore, probably inadequate to describe the conditions observed 
in multistage axial compressors. On the other hand, they provide at least a 
_ qualitative explanation of the phenomenon, and they are in agreement in their 

conclusion that steadily rotating patterns may occur in stalled, rigid blade rows. 

In addition to the limitations imposed by their simplifying assumptions, 
_ these theories all have the drawback that they involve several unknown para- 
meters. In the analysis based on ‘airfoil’ relations, for example, both the lift- 
curve slope m and the lag angle A are involved. Unfortunately, one’s under- 
standing of stalled-airfoil flow is insufficient to determine either of these 
quantities with assurance. Each of the analyses discussed has this drawback, 
although the governing parameters are expressed in other terms. It appears 
that more definitive theoretical and experimental work on stalled airfoils and 
(or) channels is needed. 

Presumably, the several different analyses, being alike in basic structure but 
differing in the assumed blade-characteristic relations, have relative advantages 
and disadvantages, but it is difficult to choose between them on the basis of 
the available experimental data. The present author believes that there is a 
considerable virtue in the assumption of MENDELSON’s type of lag, i.e., phase 
lag independent of frequency, since this avoids the limitation of flow patterns 
to purely sinusoidal ones. This limitation has occurred, probably unintention- 
ally, in the work of Emmons [1] and STENNING [9]. Undoubtedly MARBLE’s use 
of a discontinuous characteristic [8] provides information unavailable in the 
linearized treatments. It seems questionable whether the blade-characteristic 
function of Figure 9 is realistic, however, for the static-pressure rise through a 
blade row surely does not vanish when stalling occurs. In this connection, the 
contrast between MARBLE’s assumptions and those of EMMons and STENNING 
is striking: the former assumes a violent reduction of static pressure down- 
stream of the stalled portions of the blade row, while the latter assume that the 
same quantity is uniformly distributed. It should not be difficult to make a 
choice between these two assumptions on the basis of suitable future experimental 


observations. 
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Addendum Supplied October 1955 


Since the manuscript of the present paper was submitted for publication, 


additional investigations of rotating stall, including interesting experimental © 
results, have appeared (References [12], [13], [14], [15], [16]). The newer ex- 
perimental results (References [13] and [16]) are particularly pertinent to any 


assessment of the theories, because they have been obtained in single-row | 
configurations, which more closely resemble the theoretical models. Also, since _ 


static-pressure fluctuations were measured beforé and behind the rotors in 


these tests, these results cast considerable light on the assumptions made by — 
MARBLE and STENNING in their theories. 


[9 


| 
. 
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The agreement between measured rates of stall propagation and those _ 
predicted by the theories of Reference [7] is similar to that shown here in 
Figure 10. (The present author does not agree with the way this comparison | 
was carried out in Reference [16] and has made his own comparison, to which — 
the preceding statement applies.) Fortunately, the question of whether rotor 
or stator was stalled, which complicated the matter in Figure 10, does not 
occur in these single-row comparisons. The authors of Reference [13] conclude, 
in part: “The measured stall-propagation rates can be made to agree with — 
theoretical predictions by adjusting certain theoretically undetermined factors 
appearing in the theories.’ 


ty 
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Abstract 


The phenomenon known as ‘rotating stall’ is described. Basically, it involves 
a nonuniform pattern of flow, steadily rotating relative to both the fixed and the 
- rotating blades of axial-flow compressors. Attempts to analyze the phenomenon 
_ by means of small-perturbation theories are reviewed. It is shown that the work 
of several investigators can be included in a single formulation and that they 
' differ only in the type of blade-characteristic relations assumed. The present 
- author’s analyses based on linearized ‘airfoil’ and ‘channel’ relations are pre- 
sented in detail. The various theories coincide in the qualitative conclusion that 
self-induced rotating-stall patterns can occur in stalled, rigid blade rows. A 
: comparison with available experimental results is inconclusive, so far as a choice 


between the several theories is concerned. 


OR eT Se ee LD 


z Zusammenfassung 


Die als «umlaufende AbreiBstromung » bekannte Erscheinung wird beschrie- 
ben, welche durch ein Strémungsbild gekennzeichnet ist, das sowohl beziiglich 
des Stators wie auch des Rotors in einem axialen Verdichter stetig umlauft. 
Fine Ubersicht der theoretischen Ansatze zur Beschreibung dieser Erscheinung 
im Rahmen einer Stérungsrechnung ist gegeben. Es kann gezeigt werden, dass 
nsdtze einheitlich behandelt werden konnen und sich nur in der 
angenommenen Schaufelcharakteristik unterscheiden. Berechnungen des Verfas- 
-sers auf Grund von linearisierten «Auftriebskurven» und « Kanalbeziehungen » 
sind ausfiihrlich dargestellt. Alle Theorien liefern das gleiche qualitative Ergebnis 
eines selbstinduzierten, umlaufenden Strémungsbildes in einem starren Gitter 
mit abgerissener Stromung. Ein Vergleich mit Versuchen ergibt vorlaufig keine 
zwingenden Schliisse betreffend die beste Wahl zwischen den Theorien. 
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Uber den Fehler des Runge-Kutta-Verfahrens fiir die 
numerische Integration gewohnlicher Differentialgleichungen 
n-ter Ordnung 


Von WALTER GAUTSCHI, Basel!) 


Die klassischen Runge-Kutta-Formeln zur numerischen Integration eines 
Systems von gewéhnlichen Differentialgleichungen erster Ordnung kénnen im 
Prinzip auch fiir Differentialgleichungen héherer Ordnung benutzt werden, 
indem eine Differentialgleichung n-ter Ordnung in bekannter Weise als ein 
spezielles System von n Differentialgleichungen erster Ordnung aufgefasst 
wird. Den entsprechenden Approximationsfehler hat BrEBERBACH [2]?) mit 
Hilfe der Taylor-Entwicklung abgeschatzt®). Indessen sind in neuerer Zeit von 
ZURMUHL [7], [8]*) allgemeinere Runge-Kutta-Formeln aufgestellt worden, die 
die numerische Integration von Differentialgleichungen n-ter Ordnung auf 
direktem Wege durchfiihren. Die Genauigkeit der Approximation erfahrt 
dadurch, besonders fiir 7 > 2, eine bemerkenswerte Steigerung [7]. Im fol- 
genden soll auch zur Zurmiihlschen Fassung der Runge-Kutta-Formeln eine 
Fehlerabschatzung (fiir den ersten Integrationsschritt) mitgeteilt werden®). 


Die Runge-Kutta-Formeln fiir Differentialgleichungen n-ter Ordnung 


Sei die Differentialgleichung in der nach der héchsten Ableitung aufgelésten 


Form 
Mee Nea pa) (1) 


gegeben und an der Stelle « = x) das System von Anfangswerten 


a) (264) ar ya (v == 0, ce has oa Mie 1) (2) 


Mathematische Anstalt der Universitit. 
Die Zahlen in eckigen Klammern verweisen auf das Literaturverzeichnis S. 461. 


3) Von M. Lorxin [5] ist eine weitere, aber unvollstandige Fehlerabschatzung durchgefiihrt 
worden. 


4) Vel. auch [4], S. 28. 

5) Inzwischen sind fiir 2 = 1 von ALBRECHT [1] und fiir beliebiges » von BuKovics [3], aus- 
gehend von anderen Voraussetzungen, weitere Fehlerabschatzungen veréffentlicht worden, die im 
Prinzip auch die Fehlerfortpflanzung beriicksichtigen. Uber den Einfluss der Fehlerfortpflanzung 
vgl. auch eine Bemerkung von RurisHAUSER [6]. : 


) 
-) 


ee 
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_ vorgeschrieben. Wir setzen!) 


_ h=% = Xo, 


% Te Sarete fyn-?-1 
F(a) = 96) + ft P+ +e ar? =O. .m- 1), 
Ry Noigy 1 ylO)s = sk 40), «2-1, ON | 


: bao froe 2 Ba) + aan bo Ta) 


4 peat pve = T,_s(>) ue hy], 
b= trots. Ts) t+ gee te Ta) (3) 


An-e# 


hn 
i [xo +h, Ty) + — hy, F,(2) 
VMs es 
7 oT eae + hhs| | 
und ferner 
, 4 
b= DS yyoh ¥=9,1,...,0-1) (4) 
e=1 
mit 
Pei Sam Boo aed SD ae 
ce ee oe) ee eA ce 


ee 6) 
Ve eae Ps | 

Durch das Runge-Kutta-Verfahren werden nun fiir die »-ten Ableitungen 
y)(x) der gesuchten Lésungsfunktion in der folgenden Weise Naherungs- 
funktionen y°)(x) definiert : 


y)(4_ + h) = T(1) + ST Pers (lens, eo (6) 
Voraussetzungen 
Wir definieren » Funktionen /, der 1 + 1 Variablen x,-7,,.--, Yn durch 


h= us (A= 1, Bytes. 15 fa=T(%, Mas ces at 


1) Die folgenden Bezeichnungen sind gegeniiber den Darstellungen in [4] und [7] leicht geandert. 
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Es sei R ein Bereich 
(R) le *%| Sa, (4, - "184 C=1,2, 22,9) 


des (x, ¥4,.-., V,)-Raumes, in dem die Funktion f(x, Vy» +++, Yn) und. ihreg 
Ableitungen 0'f/0x' dy --- dyin ,0 <1 S< 4, eindeutig und stetig sind. Wir set- 
zen voraus, dass in R 


[ye NE ial 2 ae (7), 
| i ee M J 
/—— —| <= ———____ < 8 
| ox" ee = Nit tin-t (O<1 54) ( 1 
t 
gilt!) und dass ferner die Ungleichungen ; 
| 
an? av n ad 5 
Not = on (ea OT ena, ed), eee at = (9)" 


erfiillt sind. Die ersten » Ungleichungen in (9) sorgen unter anderem dafiir, 
dass fiir |* — x,| <a die Werte [x, y(x),..., y"~"(x)] der exakten Lésungs-_ 
funktion und ihrer Ableitungen innerhalb des Bereiches R bleiben. Ausserdem 
haben sie, wegen (7), zur Folge, dass |k,| < N (1 <o <4) ist. Die letzte 
Ungleichung in (9) erweist sich im Laufe der Fehlerabschatzung als zweck- 
massig. 


Abschatzung des Fehlers 
Sel 
é,(h) = ye (4x9 + h) — y(xg +h) (v=0,1,...,"%— 1) 


der Approximationsfehler der y-ten Ableitung. Aus der Taylor-Entwicklung 


’ hn-» An-v+. af | 
(v) = as ie _ 5 eek is 
YO (to + A= Tl) + (n — »)! treo + (71—v+1)! dh ipo 


und aus (6) ergibt sich 


An-v An-» An-vr.1 aft | 
> h = — k —_ ——— — es = e i 
€,(") (n — v)! y (mn — v)| ‘h=0 (m—v-+1)! dh\rq-o 


Es ist bekannt [7], dass fir m= 0,1,...,n—»+2undn Zane 


qm 


dhm €,(h) Per ay (v = 0, Boxes ph 2) 


und fir m= 0,1,...,4 und 21 


qm 
dh™ En _(h) eet = 0 

1) (7), (8) sind im wesentlichen die Bieberbachschen Voraussetzungen, bis auf den Unterschied 
dass hier die Abschatzung (8) nur fiir die Funktion fy, und nicht fiir samtliche fa gelten muss. 


? 
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‘ist. Daher folgt aus (10) 


> aie 


kT SS ei at Fee OL ee 
Te cere as | ee 
. (222,05 yS2-2), 
a 
: frie, oth, es ge 

el = (40 ae See as ie = 1) ae 


_ fiir gewisse #; mit O< 8; <1. Aus (11), (12), (4) und (5) erhalt man fiir » = 2 
die Abschatzungen 
jh [rts 


[e,(2)| S 6 (p+ 3)! 


3 | + max|,’|) 


+ p (6 max| hy’ | + 6 max| ky’ | + max| k7’|) 
— 6 (max| &/’| — max|7”|) } 
(p= nm— wv; v= 0, by... ,% = 2) 
cond fiir 1 = 1 die Abschatzung 


h)| S a 7 {10 max| &IY| + 10 max| Al’| + 5 max|hj"| + 6 max|/""]}. 


(14) 


; le, a 


Striche bedeuten hier Ableitungen nach h. Ferner sind alle Maxima fiir |h| Sa 


- zu nehmen. 

Fiihrt man die Differentiationen in (13) und (14) aus, und schatzt man die 
- gewonnenen Ausdriicke mit Hilfe von (7), (8) und (9) ab, so erhalt man nach 
-einiger Rechnung folgende Resultate: Es gilt fiirn = 2 . 


Ol S spar OPE) + PH lm) — Kem) Ar? | 


é, 
(pb=n—v;v=0,1,...,%— 2) | 
mut 
G(n) = 2, {737 n® + 3555 n® + 3671 n + 361 + 144 44(n) (w+ 1)} MN 
ee (273 n2 + 642. + 297) M@N+ > (n+ 1) MPN, 


H(n) = 2 {965 n® + 4599 2 + 4435 m + 69 + 144 da(n) (1 + 1) } MN 


| 
| 
/ (16) 
| 
| 


ot a (369 22 + 858 n + 369) M?.N + : (n+ 1) MPN, 
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K(n) = <7 { 3354 n® + 15804 n® + 18198 n + 4578 + 864 d(x) (x + 1) }M N- 


+ se (966 n® + 2244 m + 1182) M2 N +3 (n+ 1) M?N 
und fiirn =1 . 


184320 
— 11107 + 32 6,(n) (2050 n* + 8022 n2 + 6093 n + 31) 
+ 24 6,(n) (1060 2 + 2379 m + 523) + 3840 65(n) (n + 1)} : 


M? N 
11520 


len (A) Ss Feen {42401 n* + 237440 n® + 336570 n? + 73112 
: 


a 


{7721 n> + 27621 n? + 23163 n + 3263 t (17) ; 


+ 6 64(m) (805 n? + 1712 n + 667) + 576 9(n) (n + 1) } 


MN 
960 
M+*N 
+ Taq +] LAS, 


wenn allgemein 


a ~ {293 n? + 604 2 + 275 + 36 6,(n) (n + 1)} 


— 


0 (n<k) 
0,(”) = (18) 
1 (n> k) 


gesetat wird. 
Fiir n = 1 ergibt sich aus (17) die Schranke von BIEBERBACH?) 


|v, — ¥| S[M N 3,680642361 + M2? N 5,3618055... + M3 N 1,220833... 


ses eae 


+ M*N 0,01666...]|A|> 55,37 MN |h|8 (n= 1). 


Fiir den praktisch wichtigen Fall » = 2 erhalt man 
lv. —- »| S [MN 1,086 + M2-N 0,7188 + M* N 0,025] | A|® 


<109MN To | 


is, 
| — y'| S [MN 32,57 + M? N 23,07 + M3 N 2,879 + M4 N 0,025] |h|5 


<326MN 41 |als 


(ies: 2)" 


D2 NS. 284. 
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a Man verifiziert leicht, dass fiir n = 2 

: M3—1 Ms —1 
, G(1) = ey (n+1)?>MN Moi H(n) S 20,97 (n + ON 


_ gilt, so dass man aus (15) die schwachere, aber einfachere Abschatzung erhalt: 


|e,()| S Pa (2,78 d + 3,50) (n+ 1)? MN A [alr | 


(p + 3)! it 
. (pan—9v=0,1,...,m-— 25" 22). 
f 
Ahnlich kann (17) zu 
|G,21(%)| S 0,453 (w+ 14M NA |a[s (nw 21) (17) 


vereinfacht werden. 
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Summary 


For the numerical solution of an ordinary differential equation of the m-th 
order (1) with the initial conditions (2) ZuRMUHL has developed general Runge- 
Kutta formulas which integrate this equation directly without separating it into 
a system of ” equations of the first order. In the present paper some estimations 
of the errors of this Runge-Kutta method are given, which are valid for any value 
of u. It is assumed that the function f(x, y, y’, .-.) on the right side of the differ- 
ential equation and its partial derivatives up to the fourth order are one-valued 
and continuous in.a certain neighbourhood of the initial values and that bounds 


of these functions are known. 


(Eingegangen: 7. Juli 1954.) 
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Eine genaherte Behandlung des schweren symmetrischen 
Kreisels in nicht-Eulerschen Koordinaten 


Von ERICH SPONDER, Chromepet, South India!) 


Einleitung 


; 
; 
7 


i 
4 


In zahlreichen technischen Anwendungen bildet der Kreisel ein ausschlag- 


gebendes Glied in einem mehrlaufigen Kérperverband, dessen (dynamische) 


Stabilitat hinsichtlich einer bestimmten Gleichgewichtslage sehr oft von beson- — 


derem Interesse ist. Zu deren Feststellung bedient man sich gerne der Methode 
der kleinen Schwingungen, die nur kleine Abweichungen der Glieder des « ge- 


stérten» Korperverbandes von seiner Gleichgewichtslage in Betracht zieht. — 


Unter geeigneten Bedingungen fiihrt dies zu Differentialgleichungen, die linear 
in diesen Abweichungen und deren zeitlichen Ableitungen sind, also leicht 
gelést werden kénnen, und die den Bewegungsverlauf in zwar gendherter, aber 
sehr tibersichtlicher Weise wiedergeben. Was den Kreisel betrifft, so kommt es 
dabei nur auf die Lagenanderung seiner Figurenachse an, nicht auch auf den 
Drehwinkel des ohnehin symmetrischen Rotors, und dazu gentigt schon die 
Angabe von zwei Koordinaten an Stelle von dreien des allgemeinen Falles. 

Die Methode der kleinen Schwingungen ist nur dann erfolgreich anwendbar, 
wenn die zur Beschreibung der Bewegung eingefiihrten Koordinaten zweck- 
entsprechend gewahlt werden; so versagt zum Beispiel dieser Weg, wenn die 
Lage der Figurenachse des Kreisels in den bekannten zwei Eulerschen Winkeln 
# und wy angegeben wird, wie es die Figur 1 zeigt, und es empfiehlt sich daher, 
mit anderen, besser dafiir geeigneten Winkeln zu rechnen. Ein Versuch mit 
a und # an Stelle von # und y erweist sich als brauchbar: diese Wahl ermoglicht 
tatsachlich die sehr bequeme Behandlung des schweren symmetrischen Kreisels 
mit Hilfe der Methode der kleinen Schwingungen. 

Dieses Verfahren ist nun nicht neu, die meisten Darstellungen des ange- 
wandten Kreisels bestatigen es. Doch beschrankt man sich dann nur auf den 
schnellen Kreisel und verwendet demgemiss von vorneherein vereinfachte 
Formeln. Hier soll hingegen eine anscheinend noch fehlende strenge Herleitung 
der vollstandigen Kreiselgleichungen in diesen nicht-Eulerschen Koordinaten 
gegeben werden, die dann auch fiir den langsamen Kreisel und den Grenzfall 
des spharischen Pendels giiltig sind. Dariiber hinaus ist es sogar méglich, ohne 
allzu grosse Umstindlichkeiten den Einfluss des Dampfungswiderstandes, den 


1) Madras Institute of Technology. 
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der Kreisel auf seiner Bahn um die Lotachse herum erfahrt, auf die Stabilitat 
der Bewegung zu klaren. 


Aufstellung der Kreiselgleichungen nach Lagrange 

' Die Bewegungsgleichungen des schweren symmetrischen Kreisels — gleich- 

-giiltig, in welchen Koordinaten — werden am einfachsten mit Hilfe der Lagran- 

-geschen Gleichungen zweiter Art abgeleitet. Man erspart sich dabei alle Be- 

trachtungen tiber den Drall und gelangt auf einem verhaltnism4ssig ungefahr- 
'deten Weg rasch zum Ziel. ; 

Zunachst muss die Lage des Kreisels im Raum geometrisch festgelegt 


werden. 
Lotachse 


x 
b /y] 


t Figurenachse x 
PTF 9g ; 
‘dimmer >O 
peal 


Se 2 
leltes) fy, Knotenlinie KL(Xx) 


Figur 1 
Die verwendeten Achsensysteme; alle eingezeichneten Winkel und Drehsinne sind positiv. 


Die Figur 1 zeigt das Hauptachsensystem * y z des Kreisels, der in O unter- 
stiitzt ist, und das im Raum ruhende Bezugssystem X YZ. Beide sind Rechts- 
systeme, und dementsprechend sind die positiven Drehsinne fiir die Winkel- 
geschwindigkeiten & usw. eingezeichnet. Es darf ohne Einschrankung der All- 
gemeingiiltigkeit angenommen werden, dass auch die Rotation p des Kreisels 
immer positiv sei. Die Ausrichtung des kreiselfesten Systems x y z gegentiber 
dem raumfesten System XYZ geschieht durch die drei Winkel «, und &. « ist 


f 
an Rl 
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der Winkelabstand, den die Figurenachse x von der (X Y)-Ebene hat; die Dreh- 
achse fiir « ist durch [«] angedeutet, was sinngemass auch fiir alle iibrigen’ 
Winkel gilt. 8 wird in der (X Y)-Ebene gemessen und von der Lotachse X aus” 
gezahlt. « und 6 bilden mit dem stets positiven Eulerschen Winkel ? zusammen 
ein rechtwinkliges Kugeldreieck. Der dritte Winkel, &, wird in der raumlich 
beweglichen Querebene gemessen, die senkrecht zur Figurenachse x des Kreisels 
steht, und gibt den Weg der Querachse y von der (X Y)-Ebene an. — Die be- 
sprochenen Winkel treten auch noch an anderen Stellen der Figur 1 auf und 
sind dort eingezeichnet, was zur Verdeutlichung der Darstellung dienen mag. 
— Mit #, g und x sind wie tiblich die drei rechtwinkligen Drehungskomponenten 
des Kreisels bezeichnet, und M(#) sei der Ausdruck fiir das Schweremoment, 
das um die Knotenlinie K wirkt. Dabei braucht M(@) nicht unbedingt von der 
Schwerkraft herzuriihren, es kann ganz beliebigen Ursprungs sein. Mit dieser 
Verallgemeinerung darf man also auch dann von einem «schweren» Kreisel 
reden, wenn er in seinem Schwerpunkt unterstiitzt ist und seine Gleichgewichts- 
lage nicht lotrecht ist, wie es hier in tiblicher Weise angenommen sei. — J 7 ist 
das Tragheitsmoment um die Figurenachse x, und J, ein beliebiges Quertrag- 
heitsmoment. 

Zur Anwendung der Lagrangeschen Gleichungen ist die bekannte, hier fiir 
die Koordinate « angeschriebene Beziehung 
OL OL 
Ox , O%- (1) 
fur alle Freiheitsgrade der Bewegung, also auch fiir 6 und &, auszuwerten. 
Mierin ist F, das auf den Winkel « bezogene dussere Moment, und L die gesamte 
kinetische Energie des sich bewegenden Kreisels. 

Es sei zundchst angenommen, dass das aussere Moment nur von der Schwer- 
kraft herriihrt. Da diese ein Potential besitzt, lassen sich die Momente ey 
und J; leicht als dessen partielle Ableitungen nach den Koordinaten a, Bund & 
bestimmen. Ist gm das Gewicht des Kreisels und greift es im positiv gerech- 
neten Abstand s vom Unterstiitzungspunkt O an, so erhalt man nach der 
Figur 1 fiir das Potential V, wenn es auf das durch 3 — 0 gegebene Niveau 
bezogen wird, 


era 


V=+Fems(1—cosd). 


Hiebei — und auch weiterhin — gilt das obere Vorzeichen fiir den Spielkreisel 
(Schwerpunkt tiber dem Unterstiitzungspunkt) und das untere fiir den Pendel- 
kreisel (Schwerpunkt unter dem Unterstiitzungspunkt). Mit cos? = cosa cos B 
ergibt sich dann fiir die linken Seiten der Lagrangeschen Gleichungen (1) 


OV OV 


F= Pry es SiN ReOS Bie Tigesie oF +g ms cosa sin , 
OV 


iS 
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_Die kinetische Energie L des Kreisels wird am einfachsten aus den Trag- 
_heitsmomenten und den Drehungskomponenten um die drei Hauptachsen 
_ berechnet: 


I af 
Ge ee TP) (2) 


Dazu miissen vortibergehend die #,q,7 durch a, 6, € und deren zeitliche 
Ableitungen «, 6, € ausgedriickt werden: i 


Pp=—fsina+é, qg=acosE+ fcosasingé, r=—asiné+ Bcosacosé. 


Nunmehr kénnen die Lagrangeschen Gleichungen entsprechend der Vor- 
schrift (1) zusammengestellt werden: 


I, % = +gmssinacosf — I, p B cosa — I; B* sin« cos a, 
Tragheitsmoment Schweremoment Kreiselmoment Zentrifugalmoment 
I,Bcosn =+tgmssinB +I,pa +2i,afBsin«, (3) 
Tragheitsmoment Schweremoment Kreiselmoment Coriolis-Moment 


I, p = 0 oder # = constant (aber — + constant!) 


oberes Vorzeichen: Spielkreisel; unteres Vorzeichen: Pendelkreisel. 


Damit sind also die strengen Bewegungsgleichungen fiir den schweren 
symmetrischen Kreisel in den nicht-Eulerschen Koordinaten der Figur 1 ge- 
wonnen; sie lassen sich als Momentengleichungen um die drei zueinander senk- 
rechten Achsen [a], ¢ und x deuten. Es mag vielleicht noch erwahnt werden, 
dass sich der Wirkungssinn der Kreiselmomente nach der Regel vom gleich- 
stimmigen Parallelismus bestimmt, das heisst, dass die Figurenachse x auf 
kiirzestem Wege zur Achse der Zwangsdrehung B bzw. « hin strebt. 

Wird nur eine Stérung des Kreisels aus seiner Gleichgewichtslage [woftr 

_ M(#) = 0] in Betracht gezogen, so geniigt es, die nun folgenden Untersuchungen 
auf kleine Neigungswinkel # zu beschranken. Dann sind aber auch die Winkel « 
und f klein, und es darf fiir sie angenahert der Sinus durch den Winkel und der 
Kosinus durch 1 ersetzt werden. Weiterhin ergibt sich aus den Gleichungen (3), 
dass das Zentrifugalmoment und das Coriolis-Moment nur von dritter Ordnung 
klein gegeniiber den anderen Gliedern sind, so dass beide Momente vernach- 
lassigt werden kénnen. — Ferner wird noch, um das umstandliche + zu ver- 
meiden, statt g ms der Momentenanstieg 0M/09 mit dem Neigungswinkel ein- 

gefiihrt ; je nachdem, ob es sich um einen Spiel- oder Pendelkreisel handelt, ist 
0M 0% positiv oder negativ. - Damit ergibt sich 


i OM ; - OM : 
Le =e a—-I,pp, 1,B= op Pt hPa, I,p=constant. (4) 
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Dies sind die fiir kleine Stérungen « und f giiltigen genaherten Differential-_ 
gleichungen des schweren symmetrischen Kreisels in den nicht-Eulerschen Koor- 
dinaten der Figur 1. Sie sind linear in den beiden Unbekannten « und B und 
deren zeitlichen Ableitungen, lassen sich also leicht lésen, wie es ja als Ziel der 
Anwendung der Methode der kleinen Schwingungen erstrebt wurde. 

| ; 
Lésung der Kreiselgleichungen mit einer komplexen Veradnderlichen — 

Ehe auf den Kreisel eingegangen wird, soll noch kurz der Sonderfall des 
spharischen Pendels gestreift werden. Fiir dieses ist ja die Rotation p = 0, und 
aus den Gleichungen (4) wird einfach 


# OM 


Das sind aber zwei gewohnliche Schwingungsgleichungen, die sich einzeln 
lésen lassen. Mit 0/09 < 0 ergeben sich zwei harmonische Schwingungen mit 
derselben Pendelfrequenz 


es OM 
aarre 


fie [rad/s] . 

Bei den hier vorausgesetzten kleinen Ausschlagen kann man nun die Winkel 
a und #8 genau genug wie kartesische Koordinaten in einer Ebene betrachten, 
die sich senkrecht zur Lotachse X im Abstand 1 vom Unterstiitzungspunkt O © 
befindet (Figur 1). Aus der Zusammensetzung der beiden Pendelschwingungen 
folgt dann eine stillstehende elliptische Bahn der Pendelspitze (elliptisch-har- 
monische Bewegung). Das steht aber im Widerspruch zu ihrer wirklichen Be- 
wegung, bei der der Scheitel der etwa elliptischen Bahn standig langsam 
vorrtickt (in der Literatur oft « Prazession» genannt). Wie sich jedoch zeigen 
lasst, erfolgt dieses Vorriicken nur verhialtlich dem Quadrat des mittleren 
Pendelausschlages, so dass es natiirlich von einer Naherungsmethode nicht 
erfasst werden kann, die lediglich kleine Gréssen erster Ordnung beriicksichtigt. 
— Soviel iiber das sphirische Pendel; alle weiteren Betrachtungen gehéren nun 
wieder dem Kreisel mit von Null verschiedener (positiver) Rotation #. 

Es wird versucht, die auch sonst oft bewahrte massfreie Darstellung bei der 
Auswertung der Kreiselgleichungen (4) einzufithren. Dazu wird statt der Zeit 
t [s] eine massfreie Zeit t = ¢/T verwendet, die mit Hilfe der Zeiteinheit 


I, 
I, p [s] 


gebildet wird. Die Ableitungen des Winkels « nach der Zeit t, die immer durch 
dariibergestellte Punkte bezeichnet werden, nehmen dann folgende Formen an: 


/ " 


lod sf 
und «= 
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dies gilt ebenso fiir 6, wobei die Striche jetzt die Ableitungen der Winkel nach 


_ der massfreien Zeit t bedeuten. — Aus den Gleichungen (4) wird dann 


I, 0M/00 

[, 1 p® 
Hierin ist c der einzige noch verbleibende massfreie Parameter. Neben der 

baulich festliegenden Grésse J,/J,; enthalt er im wesentlichen den Momenten- 


ae +cat p=0, BY’ +cP—ao' =0, mit c¢ (5) 


_anstieg 0M/00 im Verhaltnis zur (doppelten) kinetischen Energie I, p? der 


Drehung um die Figurenachse, die ja den Hauptanteil der gesamten kinetischen 
Energie des Kreisels ausmacht [Formel (2)]. Das Vorzeichen von c ist dabei so 
gewahlt, dass es positiv ausfallt, wenn das Schweremoment riickfithrend ist 
(0M /0o < 0, Pendelkreisel), und negativ, wenn es umkippend wirkt (0M /00 >0, 
Spielkreisel). Ist das Schweremoment nur gering im Vergleich zur kinetischen 
Energie des Kreisels, so wird er schnell (oder stark) genannt, und ¢ ist dann 
seinem absoluten Betrag nach klein; andernfalls handelt es sich um einen 
langsamen (oder schwachen) Kreisel, wobei es wiederum nicht im einzelnen 
darauf ankommt, wie gross das Moment oder die Rotation ist. — Die Grdsse c 
charakterisiert also ihrem Betrag und Vorzeichen nach sehr gut die Belastung 


des schweren symmetrischen Kreisels und wird deshalb als relatives Schwere- 


moment bezeichnet. 

Die Glieder c « und c f nehmen in den Kreiselgleichungen (5) die Platze 
ein, wo auch sonst bei Schwingungsgleichungen die Fesselung an eine Gleich- 
gewichtslage zum Ausdruck gebracht wird, wahrend die ebenfalls massfreien 
Kreiselmomente +’ und —«’ in der charakteristischen Anordnung gyrosko- 
pischer Terme auftreten; dies aber weist auf die Méglichkeit hin, die beiden 
(kartesischen) Koordinaten « und f komplex zusammenzufassen, und davon 
soll nun hier Gebrauch gemacht werden. 

Blickt man in der Figur 1 gegen die Lotachse X zum Unterstiitzungspunkt O 
hin, und sieht man das kleine Kugeldreieck « 6 # genau genug als eben an, so 


_ bietet sich etwa folgendes Bild: 


Lotachse X 


Figurena chsex 


e 


Figur 2 
Die komplexe Veranderliche 0. 
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Die Lage der Figurenachse x wird jetzt statt durch « und # nur durch die 


einzige komplexe Koordinate F=a+i 6 dargestellt, die man ebensogut als 
einen Vektor auffassen kann. Um die komplexe Zusammenfassung selbst durch- 
zufiihren, wird die zweite der Gleichungen (5) mit 7 multipliziert und zur ersten 
hinzugezahlt; man erhalt dann 

ake aa te Oe ae. a8. 
oO" + c9—-—1%H =0 mit Oh ge und Divae ag (6) 
Wird zur Lésung von (6) fiir die Unbekannte & der Ansatz # = C exp ZT ein- 
gefiihrt, so ergibt sich zunachst die charakteristische Gleichung 


2+¢—42=0, 


und daraus folgen schliesslich die beiden Wurzeln 


i 
$ 
i 
; 


ie 
42=> 4 gic 93 (7) 


Solange c = —1/4 bleibt, bedeuten diese beiden rein imagindren Werte — 


von z zwei ungedampfte Schwingungen, die nichts anderes sein kénnen als die — 
Nutation und die Prazession. Werden ihre massfreien Frequenzen mit VV und — 


Q bezeichnet, das heisst z, = 7 N und z, = 7 Q gesetzt, so ergibt sich zu ihrer 
Bestimmung die Frequenzformel 


kaa alas 1 | 
N Gorge |p thos. 2 =p: (8) 


Die auf die Sekunde bezogene Winkelgeschwindigkeit der Nutation » [rad/s] 
und der Prazession @ [rad/s] berechnet sich einfach nach Teilung durch die 
Zeiteinheit T = I,/(I,p) zu 


pen eR eee 5 Raa 1 
none MO (Ett +4); aE eer a (9) 


Die Form 0 = C, exp 727+ C,exp i Nt (C,, Cy komplex) ist jetzt der 
unmittelbare Ausdruck fiir die beiden kreisformigen Teilbewegungen der Figu- 
renachse x, die nun nicht mehr aus je zwei Teilschwingungen der «- und f-Koor- 
dinate zusammengesetzt zu werden brauchen. Die Figur 3 mag zur geometri- 
schen Veranschaulichung dienen: 

Die Prazession kann man sich durch die Drehbewegung eines starren Armes 
um die Lotachse X erzeugt denken, wahrend die Nutation aus einer ihr iiber- 


lagerten Drehbewegung eines zweiten Armes besteht. Es ist iibrigens dabei — 


auseinanderzuhalten, ob die Nutation vom Prazessionsarm aus gezahlt wird, 
wie meist in der Kreiselliteratur (relative Nutation N,), oder von einer raumlich 
festen Anfangslage aus, wie im vorliegenden Bericht (absolute Nutation NV, oder 


Nina ES coe bi 
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Nv 


: oe Lotachse X 
absolute Nutation N 


(i gelegentlich Nj genannt) 


Prazession $2 


relative Nutation N, 


Figurenachse x 


8 


Figur 3 
Prazession und Nutation. 


auch N, zur gelegentlichen besseren Unterscheidung). — Es gilt augenscheinlich 
N,=2+4+N, [massfrei] und entsprechend »,=©@ + ¥, [rad/s], (10) 


was beim Vergleich mit anderweitig gefundenen Ergebnissen wohl zu beachten 
ist. 

Der Inhalt der Formeln (8) und (10) kann sehr itbersichtlich in einem 
Diagramm dargestellt werden: 

Fiir kleines |c| (schneller Kreisel) ist genahert 


N,=1+¢e; Q=—-¢; N, ts 2c. (11) 


Die Anzahl » der Nutationen, die auf einen Prazessionsumlauf entfallen, 
ergibt sich dann fiir einen schnellen Kreisel angenahert zu 


Se colnh 
|2| 


: se iL 
0 i 
In der Literatur wird  zuweilen als Mass fiir die Schnelligkeit des Kreisels 
benutzt; dieser Hinweis mag den Zusammenhang mit dem hier eingefiihrten c 
' geben. 
Wie die Figur 4 zeigt, ist die Nutation V immer positiv, das heisst, ihr Dreh- 
sinn ist derselbe wie der der Rotation # des Kreisels. Anders bei der Prazession 
Q: je nachdem, ob es sich um einen Spiel- oder Pendelkreisel handelt, lauft 
sie mit oder entgegen . 

Wenn ein Spielkreisel allmahlich auslauft, so gibt es fiir seine Langsamkeit 
eine bemerkenswerte Grenze, namlich bei ¢ = —1/4: dort verschwindet die 
‘urspriinglich schnelle relative Nutation vollstandig, und es bleibt nur noch 
eine viel rascher gewordene glatte Prizessionsbewegung itibrig, mit der die 
Figurenachse die Lotachse umkreist. Bei weiterhin verzdgerter Rotation # ist 
das allerdings nur ein voriibergehender Zustand; denn fiir ¢ < —1/4 erhalt nach 
dem Ausdruck (7) die Wurzel z, einen positiven reellen Anteil, und das bedeutet 
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on 


= at she relative Nutationsfrequenz, 
fae bezogen auf die 
Prazessionsbewegung 
= oe. “4 absolute Nutationsfrequenz, 
f bezogen auf den Raum 
(wie.die Prazessionsbewegung) 
Pes OMee relatives Schweremoment 
=G,29 0, Pe Pe 
kraftefreier 
Kreise/ I 
t 


a (absolute) Prazessionsfrequenz 


Figur 4 


Frequenzdiagramm. 


ein standiges Anwachsen des Neigungswinkels 3. Der Spielkreisel fallt also 
dann um, er hat die Stabilitatsgrenze bei 


2 
Li aor teva 
tiberschritten. 

Beim Pendelkreisel kann die Rotation ohne weiteres bis zum Stillstand ver- 
langsamt werden. Im Grenzfall p = 0, also c + co (spharisches Pendel), wird 
das Verhialtnis »,/m = N,/Q—=—1. Die beiden Arme der Figur 3 drehen sich 
dann mit gleicher absoluter Winkelgeschwindigkeit, aber entgegengesetztem 
Drehsinn umeinander, und die entstehende Bahn der Figurenachse x um die 
Lotachse X wird eine genaue Ellipse (bekannte Konstruktion!). Das ist die 
bereits auf anderem Wege gefundene Bewegung des spharischen Pendels. 

Der folgende Hinweis mag vielleicht noch bemerkenswert sein: 


. ek 
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Aus der Frequenzformel (8) ergibt sich, wenn der Deutlichkeit halber N, 
an Stelle von N geschrieben wird, 


ss 4 4 
oder 
P 0M/08 0M\ |oM/0d| 
2 = 

Lye ear a =T?(sen =a} 3 

ee) 

woraus folgt . 
0M 2 \ 
vg =(senGe| we; || >|@o|>|ol. (12) 


Das bedeutet erstens, dass die Nutation y, und die Prazession w beim Spiel- 
kreisel (0M/09 > 0) immer gleichen Drehsinn und beim Pendelkreisel 
(0M /09 < 0) immer entgegengesetzten Drehsinn haben; zweitens, dass die 
Prazession w immer langsamer ist als die Eigenfrequenz w, des Kreisels, wenn 
man ihn ohne Rotation (p= 0) als gewdhnliches ebenes oder sphiarisches 
Pendel schwingen lasst. Dariiber hinaus zeigt diese Beziehung aber auch noch 
einen iiberraschend einfachen zahlenmassigen Zusammenhang zwischen diesen 
Gréssen: fiir ein und denselben schweren Kreisel gilt demnach ganz unabhangig 
von seiner jeweiligen Rotation #, dass das Produkt aus (absoluter) Nutations- 
frequenz y, und Prazessionsfrequenz @ immer gleichbleibend ist. Hat man also 

_ irgendwie das w ermittelt, sei es auch aus den Kreiselgleichungen in Eulerschen 
~ Koordinaten, so folgt daraus sofort das zugehérige v, und braucht nicht erst 
durch eine gesonderte Berechnung bestimmt zu werden. 

Weiter folgt aus derselben Formel 


N,+ 2 =1; 


was in gleicher Weise auf I 
Yeto= A p (13) 
t 


fiihrt. Diese Beziehung ist — im Gegensatz zur vorigen — frei vom Schwere- 
moment M(#); sie besagt, dass die algebraische (!) Summe von (absoluter) 
Nutationsfrequenz v, und Prazessionsfrequenz m unabhangig von M () immer 
gleichbleibend ist, solange sich die Rotation p nicht andert. 


Gelegentlich ist der Rotor eines Kreiselgerates nicht von einem Gehause 
umschlossen, so dass er in freier Berithrung mit der als ruhend angenommenen 
Luft bleibt. In solchen und dhnlichen Fallen entsteht bei seiner Prazessions- 
und Nutationsbewegung ein Moment aerodynamischer Natur, der Magnus- 
Effekt. Dieser dussert sich darin, dass die Figurenachse eine seitliche Ablenkung 
auf ihrer Bahn relativ zur umgebenden Luft erfahrt, und zwar in genau dem- 
selben Sinn wie unter der Einwirkung des Kreiselmomentes. Der Magnus- 
Effekt kommt also einer scheinbar erhéhten. Rotation # des Kreisels gleich. 
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Der schwere symmetrische Kreisel mit Dampfungswiderstand 


Es soll nur der einfache Fall behandelt werden, dass der Kreisel wohl auf 
seiner Bahn um die Lotachse herum einen Bewegungswiderstand erfahrt, nicht — 


aber, dass seine Rotation abgebremst wird. Die dritte der Bewegungsgleichun- | 
gen (4), die sich auf die Figurenachse bezieht, bleibt daher ungedndert, und 7 
die ersten beiden miissen durch die Einfiihrung eines Dampfungsmomentes ver- — 
vollstandigt werden. Dieses wird — wie iiblich — verhaltlich der Bewegungsge- — 
schwindigkeit der Figurenachse angenommen, so dass seine Komponenten in : 
der «- und f-Richtung (Figur 2) in der Form ' 

OM . OM ; : 

as bzw. - B 


angesetzt werden k6énnen; darin ist 0M [oo ein konstanter, negativer Wert. — 
Damit lauten dann die ersten beiden Lagrangeschen Gleichungen (4) 


- OM oM .. : A OM oM ; 3 
Wt Noirs ag & + oa a—I,pp, I,p= 8 p+ oa B+iI,pa. (14) 


Wird wiederum die bequemere massfreie Darstellung benutzt, so andert 
sich gegentiber den Gleichungen (5) nichts, als dass noch ein Glied mit «’ bzw. 
B’ hinzukommt, das mit dem Faktor 


T 0M OM /08 
= ee ee = a 15 
I, od Ip 1) 
multipliziert ist. Dieser immer positive Ausdruck a kann als ein Mass fiir das 
relative Dampfungsmoment angesehen werden und sei kurz mit Dampfung 
bezeichnet. — Damit erhalten die Gleichungen (14) des gedampften schweren 


Kreisels fiir kleine Neigungswinkel die genaherte, massfreie Gestalt 
a” +aa'+ca+t p’=0, B'+ap'+cB—«'=0. (16) 


Obwohl eine sofortige Lésung dieser simultanen, linearen Differentialglei- 
chungen ohne weiteres méglich wire, soll auch hier vor der weiteren Behand- 
lung deren vorteilhafte Zusammenfassung unter Verwendung der einzigen 


komplexen Verdnderlichen @ = « + ? B vorgenommen werden: die Bewegung 
der Figurenachse x erscheint dann in der vektorischen Form 
O +08’ +c8—i8' =0. (17) 


Wird jetzt der iibliche Lésungsansatz arn exp 27 eingefiihrt, so ergibt 
sich die komplexe charakteristische Gleichung 


2Z+az+c—iz=0 (18) 
mit den zwei Wurzeln 


ae=(-F+5)4 hike *: bt be: c. (19) 


hm} a) 


. 
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-Beide Werte sind nunmehr komplex, was besagt, dass die durch sie be- 
schriebene Nutation und Prazession gedampft oder entdampft verlaufen. Zum 


-Unterschied vom spharischen Pendel kénnen aber hier beim Kreisel selbst bei 


noch so grosser Dampfung a die Wurzeln niemals rein reell ausfallen, wie es ja 


schon aus der Gleichung (18) hervorgeht. Die beiden Teilbewegungen, die sich 


nach der Figur 3 zusammensetzen, haben daher immer die Gestalt von loga- 
rithmischen Spiralen, die in abnehmendem, aber auch zunehmendem Sinn 
durchlaufen werden kénnen. Um dies festzustellen und um die (massfreien) 


_ Frequenzen von Nutation und Prazession zu bestimmen, braucht man nur nach 


den Rechenregeln fiir komplexe Zahlen die reellen und imaginaren Anteile von 
z, und z, aufzusuchen; es ist ja allgemein 


Va+ib=V(Var+ 2+ a)i2+i V/(Vart 08 a)/2. 


Dieses Verfahren ist aber etwas schwerfallig, und so mag noch eine zeichnerische 
Lésung wegen ihrer Einfachheit angegeben werden. Sie besteht darin, dass die 
beiden komplexen Zahlen z, und z, nach der Formel (19) in der Gaufschen 
Zahlenebene geometrisch konstruiert werden. Die Figur 5 zeigt die Durchfiih- 
rung dieses Lésungsgedankens fiir einen gediémpften Spielkreisel mit den 


qz, Mutation 
\ 


Achse 


ee, Prazession 


0,2 reelle 
Achse 


(0,103) +42 


Figur 5 : 


Die Konstruktion der komplexen Wurzeln 2, und 2). 


4,2— (-01+ 5) Le /(- 0,1 + Ais aie yee 


Zunachst wird die Zahl {— 0,1 + (z/2)] eingetragen und zur Konstruktion von 
Ulm Ol oe (2/2)]? + 0,2 quadriert (verdoppeltes Argument, was durch die stark © 
eingezeichneten gleich grossen Drehwinkel um O angedeutet ist, und auf das © 
Quadrat gestreckter Betrag) ; dann wird noch + 0,2 hinzugezahlt und aus dem © 
so erhaltenen Ausdruck die Wurzel gezogen (im umgekehrten Sinn wie das t 
Quadrieren). Das Zwischenergebnis Vi= 0,1 + (2 /2)] + 0,2 wird schliesslich zur — 
Zahl [—0,1 + (z/2)] geometrisch addiert und von ihr subtrahiert, womit die — 
beiden Wurzeln z, und 2, gefunden sind (es ist z= —0,28,+ 0,77,7 und : 
2 = + 0,08, + 0,22, 2). 3 

Zur Probe kann noch die geometrische Addition von z, und z, ausgefiihrt 
werden, die ja den Wert (—a+ 7) haben soll, und ferner die Multiplikation — 
%, ,, die gerade c ergeben soll. Dies letztere erfordert, dass das Produkt der — 
beiden absoluten Betrage |z,| |z,| gleich |c| ist und dass die Summe der Argu- 
mente von z, und z, genau 2 Rechte (fiir c<0) bzw. 4 Rechte (fiir c>0) _ 
betragt. Aus dieser Bedingung ist iibrigens sofort zu schliessen, dass fiir einen — 
Spielkreisel (c << 0) unméglich beide Wurzeln in der linken Halbebene legen 
kénnen, so wie es zum Beispiel aus der Figur 5 fiir die Prazession mit ihrem 
positiven Realteil ersichtlich ist. — Eine in der rechten Halbebene liegende © 
Wurzel bedeutet aber eine Entdimpfung der Kreiselbewegung. Danach ware 
also der Spielkreisel mit Dampfungswiderstand instabil, und das ist ganz 
offenbar im Widerspruch zu der alltaglichen Erfahrung, die ja gerade im Gegen- | 
teil eine ausgesprochene Stabilitat der aufrechten Lage zeigt. Nun, diese Er- 
scheinung des Aufrichtens wird durch eine ganz andere, hier nicht beriicksich- 
tigte Ursache hervorgerufen, namlich durch die Reibung der Kreiselspitze auf 
ihrer Unterlage. Uberwiegt dieser Einfluss, so kann tatsichlich die meist nur 
schwache instabilisierende Wirkung der Wurzel 2, dem Spielkreisel nichts 
anhaben: er bleibt aufrecht. 

Beim Pendelkreisel mit c>0 liegen hingegen beide Wurzeln z, und z, 
immer in der linken Halbebene, wie auch spater noch gezeigt werden wird. 
Jedenfalls aber miissen sie so angeordnet sein, dass die Verbindungslinien O z, 
und O 2, spiegelbildlich zur reellen Achse sind; denn nur dann kann die Summe 
der Argumente von z, und z, 4 Rechte betragen. Es gilt daher einfach aus 
geometrischen Griinden fiir die Realteile (R) und Imaginarteile (J) der beiden 


Wurzeln 
|R(z,) | ey |R(zs) | 
1da)t ~ ote) (20) 
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Werten a= 0,2 und c = — 0,2, also fiir 
s 
; 


Wie die Figur 5 zeigt, besteht dieselbe Beziehung auch fiir den Spielkreisel, 
und zwar weil die Verbindungslinien O z, und O %, jetzt wegen c < 0 spiegel- 
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bildlich zur imaginaren Achse liegen miissen. Die Anteile J bedeuten aber 
nichts anderes als die (massfreien) Frequenzen von Nutation und Prazession, 
und die Anteile & sind deren Dampfungen (bzw. Entdampfung bei der Pra- 
zession des Spielkreisels). Die Beziehung (20) sagt also aus, dass die meist 
wesentlich schnellere Nutation viel starker gedampft ist als die Prazession, was 
“man ja bei jedem Kreiselversuch deutlich genug beobachten kann. Ist der 
Kreisel jedoch langsam, so behauptet sich aus eben demselben Grunde auch die 
_Nutation langere Zeit neben der Prazession, und im Grenzfall des spharischen 
Pendels sind diese beiden Teilbewegungen auch hinsichtlich des Abklingens 
‘vollkommen gleichwertig. 
Handelt es sich um einen schnellen Kreisel (kleines c), so zeigt die Kon- 
struktion der Figur 5 oder auch die Formel (19), dass sich z, auf den Ursprung O 
zuriickzieht. Bei der geometrischen Addition von z, und z, tiberwiegt dann z, 
bei weitem das wesentlich kleinere z,, und es gilt daher in erster Anndherung 


zi=(-a+1). 


Wegen 2, = c/z, ergibt sich aber daraus auch eine Naherung 


Cc Ue 
a Gata)’ 


Mit diesem Wert kann nun bs verbessert werden: 
ret . I 
z= (-at+1)—&%&, 


woraus wiederum ein noch genaueres 2}! = c/z}' folgt, und so weiter, bis mit 

"Hilfe dieses gut konvergierenden Iterationsverfahrens die gewiinschte Genauig- 
‘ keit erreicht ist. 

Im iibrigen zeigt sich fiir wirkliche Kreisel mit gewohnlich nur kleinem 

- Dampfungswiderstand, dass die beiden Wurzeln 2, und z, sehr nahe der imagi- 

naren Achse liegen. Dann ergibt sich aber die Méglichkeit, das Problem des 

‘schwach gedaimpften Kreisels genahert zu lésen, und das soll im nachsten 
-Abschnitt versucht werden. 


Eine Naherungslésung fiir kleinen Dampfungswiderstand 


Die ungedampfte Kreiselbewegung wird durch zwei rein imaginare Wurzeln 
z, und z, beschrieben, die — wie im vorigen Abschnitt ausgefiihrt — beim Hinzu- 
treten einer Dampfung a komplex werden [Formel (19)], also einen reellen Anteil 
erhalten. Ist die Dampfung a nur klein, so wird auch dieser Zuwachs nur klein 
sein, und die komplexen Wurzeln lassen sich dann bei Beschrankung auf das 
erste Taylor-Glied durch 
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ausdriicken. Hierin werden die Wurzeln z als Funktionen der Dampfung a an- 
gesehen, und (dz/da),_) bedeutet den (im allgemeinen komplexen) Differential 
quotienten an der Stelle a = 0, die ja der Ausgangspunkt zur Gewinnung der 
Naherungslésung fiir kleinen Dampfungswiderstand ist. - Wird zur Ermittlung 
von dz/da die komplexe charakteristische Gleichung (18) nach a abgeleitet, so 
erhalt man zundachst i Pana 
2 Bn ter Gag itt bea 

und daraus 

dz —Zz 

ag oe re are 


0 nindmnpmenapatalidialinaia, 


An der Stelle a = 0 sind die beiden Werte z durch die Formel (7) oder auch 
(19) gegeben, und dies eingesetzt fiihrt auf die fraglichen Differentialquotienten 


1 : rh eg rie : 
(421) Meee Coa Ge eS 
da =O = ? 
: 7+ 24 aca ay vagy 2 +e 
4 + 
f (21) 
Laie rae oe ine 
(=) ig nd 4 Serle 4 
“Ga fas rpm a 


Man bemerkt zunachst eine erfreulich einfache Form dieser beiden Aus- 
driicke. Weiterhin fallt auf, dass sie rein reell sind, so dass ein kleiner Zuwachs 
(dz/da),.9 @ senkrecht zur imaginaren Achse steht (Figur 5). Das besagt aber, 
dass die imaginaren Anteile der Wurzeln z, und z, dieselben bleiben, oder mit 
andern Worten, dass die dadurch bestimmten Frequenzen WV und Q in erster 
Naherung von einer kleinen Dampfung a unabhangig sind. Diese in der Schwin- 
gungslehre sehr allgemeine Tatsache gilt also auch hier und tragt wesentlich 
zur Vereinfachung der Lésung bei. Fiir die Frequenzen von Nutation und Pra- 
zession gilt daher bei ungedampfter sowie bei schwach gedampfter Bewegung 
gleicherweise die Frequenzformel (8). Diese lasst sich ubrigens noch benutzen, 
um die Differentialquotienten (21) zu vereinfachen; es ergibt sich in leicht er- 
sichtlicher Weise 


az, 25 ’ Tt dzy i Q 
(i ).-07 WB und (2) =+a-@ 


worin auch fiir VW — 2 nach der Gleichung (10) N, geschrieben werden kann. 
Damit werden nun die reellen Anteile der beiden komplexen Wurzeln (z,) 
und (z,)_ bei kleiner Dimpfung a 
N 


Q 
eats 5) [On a bzw. + N_— oO a 


a 
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| mit 
N20, N-Q=N,>0, QS0. 

Es zeigt sich wiederum, dass der letztere Anteil fiir 2> 0, also fir die 
Prazession des gedimpften Spielkreisels, positiv ausfallt, die Bewegung der 
_Figurenachse somit instabil ist. - Zusammengefasst lauten dann endlich die 
-genaherten Wurzeln z, und z, der schwach gedaémpften Kreiselbewegung 


Q 


a+iN(Nutation), z= + N—Q 


N : : 
= —y a +1Q (Prazession). (22) 


Man findet hier iibrigens die Beziehung (20) wieder: das Verhaltnis des reel- 
len Anteils zum imaginaren ist fiir beide Schwingungen dasselbe, und zwar 
a|(N — Q) (wenn vom Vorzeichen abgesehen wird). 

So iiberraschend einfach die Formeln (22) auch sind, man sieht ihnen nicht 
ohne weiteres an, mit welcher Genauigkeit sie das Ergebnis liefern. Insbesondere 
wird man im Zweifel dariiber sein, bis zu welchem Wert der Dampfung a man 
nur gehen darf, um nicht eine zu grosse Unsicherheit dieser Naherungen in 
Kauf nehmen zu miissen. Eine Fehlerabschatzung ist also sehr erwiinscht; sie 
kann so geschehen, dass die Summe und das Produkt der beiden Wurzeln 
gebildet und mit den Sollwerten verglichen werden, welche wegen (18) gleich 
—a+ibzw.c sind: 

Die Summe ergibt 

tm t ati (Nt Q)=-ati, 
a) 


genau wie es sein soll, und das Produkt ergibt 


2, 29 -NQ (xa) 
eas c 


Hier kommt statt c ein grosserer Wert heraus; der dabei begangene relative 
Fehler ist 


a i ae ; (0M /08)2 


a eae = (23) 


~ 146 (5p)? — 41, (0M/08) * 
Dieser Ausdruck gibt zwar nicht im einzelnen dariiber Aufschluss, wie sich 
der Fehler auf die beiden Wurzeln verteilt, doch gibt er jedenfalls einen guten 
Anhalt zur Beurteilung des Ergebnisses. So ware es zum Beispiel nicht ratsam, 
sich blindlings auf ein genahertes Ergebnis zu verlassen, das mit den Werten 
a= 0,2 undc = —0,2 erhalten wird (zeichnerisches Beispiel der Figur 5); denn 
jetzt ware schon das Kriterium (23) 
a2 0,04 


id Ses ae fo) 
aera = 1 — i078 20% 
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eine deutliche Warnung zur Vorsicht. Es ware dann unerlasslich, eine Probe: 
anzustellen und gegebenenfalls eine Verbesserung vorzunehmen, wenn man 
nicht gleich von Anfang an die bequemere zeichnerische Lésung vorziehen 
sollte. : 


Summary 


The gyroscope has many technical applications as an essential member in a 
system with several degrees of freedom, the (dynamic) stability of which with 
respect to a certain equilibrium position may be of particular interest. To investi- 
gate the stability of a system by the well-known method of small oscillations, the 
movements of the different members must be described by convenient coordinates. 
As to the gyroscope, it is necessary for this purpose to replace EULER’s angles 
by other coordinates, with which the rigorous equations of its movement are 
deduced. Assuming only small displacements of the axis of symmetry from its 
undisturbed equilibrium position, and by making the problem dimensionless, 
a simple solution can be obtained for the nutation and precession of a heavy 
Syroscope both as a symmetrical top and as a gyroscopic pendulum. Further, 
the damping of the gyroscope by its movement around the equilibrium position 
has also been considered in order to be able to estimate its influence on the 
stability. 


(Eingegangen: 25. Mai 1954.) 
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Uber den Mechanismus einiger Elektrete 


Von Karri ANTENEN, Ziirich!) 


I. Einleitung 


Die Bezeichnung «Elektret» fiir ein «permanent» polarisiertes Dielektrikum 
geht auf Ecucui zuriick. Er polarisierte Scheiben von Karnaubawachs in elek- 
trischen Feldern von 1000 bis 10000 V/cm bei héherer Temperatur und fand, 
dass im ersten Falle nach Abschalten des Feldes bei Zimmertemperatur und 
Kurzschluss der Elektroden die der Anode gegentiberliegende Dielektrikumober- 
flache negative, die der Kathode anliegende Oberflache positive Ladung zeigte. 


Bei mittleren Feldern von 10000 V/cm und mehr war der Sachverhalt wahrend _ 


weniger Tage derselbe, doch kehrte sich das Vorzeichen der scheinbaren Ober- 
flachenladung spater um. Die erste Ladungsart nennt man Hetero-, die zweite 
Homocharge. 

Es wurden Elektrete aus den verschiedensten Substanzen hergestellt, wobei 
das Schwergewicht dauernd auf polare Substanzen gelegt wurde (Wachse, Harze, 
Polyvinylchlorid, Plexiglas usw.). Dies prajudiziert eine unserer Auffassung nach 
unfruchtbare Dipoltheorie. 


') Institut fiir technische Physik der ETH. 


a. 2 eee ‘ 
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4 EwincG?) berichtet iiber R6ntgenaufnahmen von Karnaubawachs-Elektreten; 

_ die orientierte Dipole zeigen sollen. Uber Raumladungen in Wachsen berichten 

_ THIESSEN, WINKEL und HERRMANN sowie JAGER. 

: Gross?) stellte Experimente tiber das Zustandekommen der Homocharge an 
und erklarte sie als Luftionen, die durch elektrischen Durchschlag auf die Dielek- 

_ trikumoberflache gelangen. Die diesbeziiglichen Experimente sind so tiberzeu- 
-gend, dass in bezug auf die Herkunft der Homocharge kaum ein Zweifel még- 
lich ist. 


Erklaévungsméglichkeiten fiiy die Heterocharge 


: Piezo- und pyroelektrische Erscheinungen, wie sie GEMANT®) in Betracht zog, 
um die Ladungsumkehr erklaren zu kénnen, wollen wir ausser acht lassen, da 
bisher an Wachsen, Harzen usw. keiner dieser Effekte mit Sicherheit nachge- 
wiesen werden konnte. 

Damit gibt es fiir die Heterocharge grundsatzlich noch drei mégliche Mecha- 
nismen, die einzeln oder in Verbindung miteinander den Effekt hervorbringen 

_ kénnten: 

a) Dipolovientierung. Man kénnte sich vorstellen, dass sehr lange Dipole, wie 
organische Molekiile sie vielleicht hervorbringen, in kurzgeschlossenem Zustande 
hohe Relaxationszeiten aufweisen. Oder man kénnte an Dipolgruppen denken, 
die in kurzgeschlossenem Zustande aus energetischen Griinden orientiert bleiben. 
«Elektrete» in diesem Sinne konnten zum Beispiel aus Bariumtitanat hergestellt 
werden. 

-b) Grenzflachenpolarisation (J.C. MAXWELL). Unstetigkeitsflachen im Dielek- 
 trikum bauen wegen Leitfahigkeitsunterschieden Flachenladungen (das heisst 
sehr diinne Raumladungsschichten) auf. Diese erzeugen ein «Dipolmoment». 
Namentlich zwei wichtige Grenzflachen sind immer vorhanden: jene zwischen 
Elektroden und Dielektrikum. 

c) Ionenwanderung iiber makroskopische Distanzen. Dadurch kénnten an den 
Elektroden Raumladungsschichten entstehen, die einen makroskopischen Dipol 
erzeugen wiirden. Es wurde jedoch mehrfach die Vermutung gedussert, dass 
damit die «Permanenz» der Ladung infolge der endlichen Leitfahigkeit des 
Dielektrikums nicht erklart werden kénne (zum Beispiel Gross). Es wurde aber 
offenbar ausser acht gelassen, dass auch hier die auftretende Relaxationszeit 
durchaus nicht (e €,)/¢ zu sein braucht, wobei o etwa die mit Wechselspannung 
gemessene Leitfahigkeit bedeuten kann. 

Es wurde in der vorliegenden Arbeit versucht, Klarheit zu schaffen liber das 
Zustandekommen der Heterocharge. 


II. Potentiale 


Wichtige Aufschliisse kann der Potentialverlauf im Innern des Dielektrikums 
ergeben. Es wurden deshalb scheibenformige Proben aus filtriertem Karnauba- 
wachs hergestellt mit einem Durchmesser von 75 mm und einer Dicke von 12 mm, 
Die parallelen Flachen wurden mit Silber bedampft und langs des Umfanges 


Sonden eingeschmolzen. 


1) M. Ewine, An X-ray Study of the Structure on Electrets, Phys. Rev. 36, 378 (1930). 

2) B. Gross, On Permanent Charges in Solid Dielectrics I, Phys. Rev. 67, 253 (1945). — On 
Permanent Charges in Solid Dielectrics II, J. chem. Phys. 17, 866 (1949). 

3) A, GEMANT, Recent Investigations on Electrets, Phil. Mag. 20, 929 (1935). 
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An die versilberten Oberflachen wurde eine Spannung von 1200 V gelegt, 
worauf mit einer Kompensationsschaltung die Potentialdifferenzen zwischen 
Sonde und Elektrode oder Sonde und Sonde gemessen wurden. Auf diese Weise 
wurde der Potentialverlauf bei 70°C, 50°C und 35°C festgestellt. Es zeigte sich 
folgendes: 

Die Feldstarke am Rande nimmt mit der Zeit zu. Es gibt ein ausgedehntes 
Mittelstiick, in dem die Feldstarke zeitlich abnimmt, radumlich aber konstant 
bleibt, das heisst in dem keine Raumladungen auftreten. Bei 70°C betragt die 


Potential { 100 V] 


Ss 


hy So oe Oy SS Sa | se 


tf PPTs AOR PT es t : ley (s iae  aa 
mm Entfernung 
von der Anode 


Figur 1 


Potentialverteilung in Karnauba bei verschiedener Leitfahigkeit: 
TO Se 0s e* OSU om te 11 Ostag = 4:10-15 Q-1 cm}; T = 70°C. 


Feldstarke am Rande des Dielektrikums 10 000 V/cm, in der Mitte im stationdren 
Zustand 330 V/cm. 


Bei niedrigerer Temperatur ist die Randfeldstarke (durch Extrapolation 
bestimmt) geringer (T = 50°C: 8000 V/cm). Dafiir erstrecken sich die Inhomo- 
genitaten etwas weiter ins Innere des Dielektrikums. 

Die Relaxationszeiten der Potentialdifferenzen betrugen am Rande (0-0,7 mm 
von den Elektroden entfernt) bei 70°C 1-1,2 min (Stationarstrom 10-® A) und 
bei 50°C 27 min (Stationarstrom 3-10-11 A), Diese Relaxationszeiten sind ver- 
gleichbar mit jenen der induzierten Ladung. Die Dielektrizitatskonstante der 
filtrierten Probe betragt 3,5. Berechnet man aus Randfeldstarke und dielek- 
trischer Verschiebung die auf einer beriihrenden Metallplatte influenzierte La- 
dung, so erhalt man 1,4 - 10-7 Clb (70°C; Elektrodendurchmesser 75 mm). 

Die Ladung kann mittels Ladungstrennung im «dissectible capacitor» nach 
Gross direkt gemessen werden (es handelt sich um die «scheinbare Oberflachen- 


ladung»). Man erhalt bei 70°C und fiir eine Probe gleicher  Leitfahigkeit 
Gr TOs Clip. 


<sccseaptiapcamattaaiaa cata iamammaaammamiias aia eal iii 
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a Fiir einen solchen Potentialverlauf k6nnen keine Dipole verantwortlich ge- 
‘: macht werden. Dagegen kann ein solcher Potentialverlauf erhalten werden, wenn 
_ man die iiblichen Gleichungen fiir lonenwanderung im elektrischen Felde fiir den 
_ stationaren Zustand lést. Sie lauten: 
: op o*p 
er — SS = cae A ts pc Seen a 
fi ot og te egne Ox eae sat) eee © 

on 0) O?n 0) 

es ey — teh ee — ag 

~ or ae le Ox ED Ne heroes 
; Wahrend des Experimentes lassen wir die Spannung konstant. Dies liefert die 


Randbedingung 


d 
fEdx=V. (2) 
0 
Es bedeuten: 
p positive Ladungsdichte, 
n negative Ladungsdichte, 
12; Feldstarke, 


b+ (b-) Beweglichkeit der positiven (negativen) Ionen, 
D+ (D-) Diffusionskoeffizient der positiven (negativen) Ionen, 
q je s und cm® neu entstehende Ladung (positiv und negativ), 
o Rekombinationskoeffizient. 
Die Gleichungen setzen voraus: 
-1. Erhaltung der Ladung, 

2. Geschwindigkeit der Ionen proportional der Feldstarke, 

3. Gaskinetische Diffusion, 

4, Massenwirkungsgesetz fiir Rekombination, hier fir eineinwertige Ionen 

formuliert. 

Im stationdren Zustand sind die linken Seiten von (1) gleich 0 zu setzen. Man 
erkennt zunachst zwei Moglichkeiten: 

1. Gleichgewicht zwischen Leitungs- und Diffusionsstrom [JAFFE})]. 

2. «Gleichgewicht » zwischen Entstehung der Ionen und Rekombination einer- 
seits, Wanderung im Felde andererseits [g — ap = 6 0/dx (E p)|[THomson, MIE, 
HERZFELD®)]. 

Beide Theorien setzen fiir die Berechnung des Feldstarkeverlaufes bestimmte 
Randbedingungen voraus und liefern in Elektrodennahe einen exponentiellen 
Abfall der Feldstarke gegen innen. Nach beiden lassen sich aus einer gemessenen 
Feldstirkeverteilung sowie aus Anfangs- und Stationarstrom die Ionenkonstanten 
bestimmen. 

Nach EINSTEIN gilt folgende Relation: 


k = Boltzmann-Konstante, T = absolute Temperatur, ¢ = Tonenladung. 
Fir T = 70°C erhalt man ~ 3- 10-2 V. Deutet man die gemessenen Poten- 
tiale nach JAFFE, so ergibt sich D/b ~ 500V, also eine gréssenordnungsmassige 


1) G. Jarré, Theorie der Leitjahigkeit polarisierbarer Medien, Ann. Ph. 16, 217 (1933). 

2) G. Mie, Der elektrische Strom in ionisierter Luft in einem ebenen Kondensator, Ann. Ph. 13, 856 
) 
) 


Karu F. HerzFe_p, Nature of Conductivity of Insulating Oils, Ph. Rev. 287 (1931). 
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Unstimmigkeit. Da die Relation auf thermodynamischen Grundlagen beruht, 
miissen wir diese Deutungsmoglichkeit verlassen, obschon JAFF& selber solche 
Unstimmigkeiten zulasst, und die zweite Méglichkeit ins Auge fassen (THOMSON, 
Mie, HERZFELD). Man erhalt auf diese Weise folgende Ionenkonstanten: 


f= 70°C: 0% = b> =6 «10 cm?/V5, 
No 10-* Clb/em* (Gleichgewichtskonzentration der Ionen), 
D 2:10-*cm?/s_ (berechnet nach EINSTEIN), 
C= 2,5 Ome BAe: 
1 = Oe OD) 23) 10=*em3s/Vs: 
: Nona Sra Oe Elb/cms, 
Dy = To 10e!¢ cme: 


ll 


Sowohl die Grésse dieser Ionenkonstanten wie der Temperaturgang sind absolut 
verstandlich. (Nach JArrEé dagegen wiirden sich zum Beispiel Diffusionskoeffi- 
zienten von 10-4 cm?/s bei 70°C ergeben.) Aus den angegebenen_ Diffusionskoef- 
fizienten wiederum lasst sich berechnen, dass der Diffusionsstrom fiir die meisten 
in Betracht fallenden Erwagungen vernachlassigt werden kann. 

Setzt man die Leitfahigkeit des Dielektrikums herab, so wird der Potential- 
verlust im Innern linearer und die Elektrodenladung geringer (Figur 1). 


III. Stréme und Ladungen 


Die Lade- und Entladestréme zeigten den bekannten zeitlichen Verlauf. 

Lade- und Entladestréme befolgen anfanglich das Ohmsche Gesetz in jedem 
Zeitpunkt, erst bei grésseren Zeiten ergeben sich Abweichungen. Es gilt iiberdies 
naherungsweise die Beziehung, dass der Anfangswert des Ladestromes = An- 
fangswert des Entladestromes + Stationarstrom (Lade- und Entladestrom posi- 
tiv gerechnet). Diese Beziehung kann theoretisch hergeleitet werden und erweist 
sich dann als giiltig, wenn die Ladungstragerzahl im Dielektrikum sich nicht 
andert. 

Es ist iiberdies zu bemerken, dass Anfangs- und Stationarwert der Leitfahig- 
keit durchaus nicht iibereinstimmen. Wir sind deshalb gezwungen, ausser den 
oben angegebenen noch beweglichere Ionen anzunehmen, deren Beitrag zu der 
Gesamtladung kleiner als 25 % sein wird, wie mittels der Temperaturabhangigkeit 
der Anfangsleitfahigkeit geschatzt werden kann. Als Ursache ftir die «perma- 
nente» Ladung kénnen sie nicht in Betracht fallen. 

Die stationire Ladung nimmt zu mit wachsender Temperatur, was sich aus 
der betrachteten Theorie dadurch ergibt, dass die Langevin-Konstante b/a, wie 
zu erwarten, mit wachsender Temperatur zunimmt. 

Tragt man die Endladestréme logarithmisch auf [log J = 7(t)], so erhalt man 
nach kurzer Zeit bei 70°C «Geraden», Man k6énnte sie so deuten, dass man an- 
nimmt, es entlade sich eine « Polarisationskapazitat» iiber den Probenwiderstand. 
Die Dicke dieser Polarisationskapazitat (die theoretisch mit der Entfernung des 
Ladungsschwerpunktes von der Elektrode zusammenfillt bei einer exponentiellen 
Feldverteilung) ergibt sich bei 70°C zu 0,25 mm. (Aus der Feldverteilung wurden 
0,3 mm erhalten.) Diese Polarisationskapazitat schwankt von Probe zu Probe 
und ist nur fiir Proben mit derselben Leitfahigkeit reproduzierbar. 

Uberhaupt muss an dieser Stelle hervorgehoben werden, dass samtliche fiir 
den Elektreten wichtigen Eigenschaften fest mit der Leitfahigkeit verkniipft 


sind. Ebenso wie fiir den Potentialverlauf gilt dies fiir die Elektrodenladung, die 
mit sinkender Leitfahigkeit abnimmt. 


ene wee 
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_ Unter der Annahme, dass die Leitfahigkeit in Elektrodennahe linear gegen die 
Elektroden hin abnimmt (in Wirklichkeit nimmt sie exponentiell ab) und im Mit- 


- telteil des Dielektrikums konstant bleibt, lasst sich bei diinnen Raumladungs- 
' schichten fiir den Entladestrom in erster Naherung eine Formel von der Form 


herleiten 
A e-i/T 
ae A Boads | 


wobei sich A und B aus den Ionenkonstanten ergeben. Es wird dabei tiberdies 
vorausgesetzt, dass sich die Leitfahigkeit am Rande anfanglich wesentlich lang- 
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Figur 2 


Entladestréme, Karnauba filtriert, Durchmesser 75 mm, Dicke 12 mm, 
T = 70°C, logJ = f(é). 


samer dndert als die Raumladung, was gezeigt werden kann. Diese Formel 
kommt den tatsdchlich gemessenen Strémen recht nahe und zeigt auch, dass die 
Betrachtungsweise mit der «Polarisationskapazitaét» fiir mittlere Zeiten ihre 
Berechtigung hat, anfanglich aber nicht. 

Vergleicht man experimentell 


t 


f Jdtund Q, 


0 


das heisst das Stromzeitintegral und die auf die Elektroden influenzierte Ladung, 
so erkennt man, dass anfanglich kein Leitungsstrom zwischen Dielektrikumober- 
flache und angrenzender Elektrode fliesst, wie zu erwarten ist. Dies zeigt die 
Existenz einer Potentialschwelle zwischen Dielektrikum und Metall auf, die aber, 
makroskopisch betrachtet, bei weitem nicht ausreichen wiirde, um etwa Tonen im 
chemischen Sinne «abzuscheiden ». Man ist deshalb zur Annahme gezwungen, dass 
Ladungstransport zwischen Dielektrikum und Elektrode nur an einzelnen aus- 
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gezeichneten Stellen mit hoher Feldstarke erfolgt [worauf auch BINDER?) hin- — 
weist]. Denkt man sich die Raumladungen durch Flachenladungen im Ladungs- 
schwer«punkt» ersetzt, so erkennt man, dass diese Flachenladungen zwar ab- — 
nehmen wegen der Leitfahigkeit des Dielektrikums, aber auch gegen die Elektro- 

den hin wandern. Die Abnahme «aussenherum» wird mit der Zeit (wegen der 
oben angefiihrten Potentialschwelle) sehr gering sein, die « Polarisationskapazitat » 

wird aber zunehmen. Fiir eine 1 yw dicke Raumladungsschicht wird folgende Ab- 

schatzung erlaubt sein: 


T 2+ RC/2 (zwei Polarisationskapazitaten in Serie) , 
Je SE MODY, : 
C/2 = e 2) F/2d = 3,5-8,9-10-14- 20-104 ~ 6-108 F, 


=O NOs = OOO Mages 


Dass nur die unmittelbar den Elektroden benachbarte Ladungsschicht fiir die 
hohen Relaxationszeiten verantwortlich gemacht werden kann, konnte experi- 
mentell gezeigt werden, indem die Potentialdifferenzen im Innern tatsachlich 
nach einigen Tagen verschwinden, wahrend die Oberflachenladung sich verhalt 
wie bereits dargelegt. 


IV. Elektrete aus Paraffin mit Beimengungen 


Es miisste also méglich sein, Elektrete (mit Heterocharge) aus einer unpolaren 
Grundsubstanz mit geringen ionalen Beimengungen zu erzeugen. Es wurde des- 
halb zu Paraffin (Smp. 90°C) Mg-Oleat beigegeben (Promille bis Prozent), und es 
zeigte sich, dass sich die so erzeugten Elektrete sowohl in Ladung wie in der 
Relaxationszeit mit Karnaubawachs vergleichen lassen. Bei einer Beimengung 
von drei Gewichtspromillen: Mg-Oleat nahm die erzeugte Ladung wahrend zwei 
Monaten nicht messbar ab. 

Wollte man die hier erzeugte Polarisation auf Dipole zuriickfiihren, so miisste 
man Mg-Oleat eine statische Dielektrizitatskonstante von mindestens 10000 zu- 
schreiben. 


Summary 


(1) It was found that the electret effect of Carnauba wax is caused by the small 
conductivity of ions in the electric field. 

(2) The constants of ions are calculated following the theory of THomson, MIE 
and HERZFELD. 


(3) Electrets were prepared with high melting paraffins and small impurities. 


(Eingegangen: 28. Mai 1955.) 


1) F. Bryper, Uber das permanent polarisierte Dielektrikum, Z. Naturforschung 6a, 714 (1951). 
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_ Mechanismus der Vereisung unterkiihlter Wassertropfen 
durch disperse Keimsubstanzen 


Von HANS RUDOLF PRUPPACHER und RAYMUND SANGER, Ziirich') 


II. Zur Deutung der Eiskeimfahigkeit der Stoffe 


A. Epitaxie als Ursache der Eiskeimfahigkeit 


In der iibersichtsweisen Erérterung der experimentellen Ergebnisse tiber die 
Eiskeimfahigkeit der Stoffe im ersten Teil dieser Notiz [1]*) ist bereits in genti- 
gender Weise auf die Schwierigkeiten hingewiesen worden, die Befunde auf Grund 

_kristalliner Verwandtschaft zu verstehen. Weder konnte eine systematische Ab- 
hangigkeit vom Gittertypus und der Kristallsymmetrie noch von den geometri- 
-schen Abmessungen auf den Spalt- und Wachstumsflachen wahrgenommen werden. 
_ Ein markantes Beispiel dafiir, dass nicht die oben genannten Beziehungen fiir die 
_Eiskeimfahigkeit im allgemeinen verantwortlich sein konnen, liefert die Verbindung 
-Kupfersulfid (CuS), die im kristallinen Aufbau keine verwandtschaftlichen Merk- 
male mit der Struktur des Eises aufweist und dennoch dieselbe Eiskeimfahigkeit 
wie Silberjodid (AglI) besitzt. 

In einer vorlaufigen Mitteilung weist R. Montmory [2] darauf hin, dass zufolge 
kristalliner Verwandtschaft Magnesiumantimon (Mg,Sb,) und Magnesiumtellu- 
rid (MgTe) bei héherer Temperatur als AgI sich eiskeimaktiv erweisen sollten. In 

der Tat findet er fiir Mg,Sb, eine Einsatztemperatur von — 4°, wahrend sie nach 
Monrmory fiir AgI unter gleichen experimentellen Umstanden leicht tiefer liegt. 
Trotz dieser Beobachtung kann aber auf Grund des gesamten bis jetzt vorliegenden 
experimentellen Materials die kristalline Epitaxie nicht herangezogen werden, 
um die Eigenschaft bestimmter Stoffe, die Vereisung unterkiihlter Wolken un- 
mittelbar einleiten zu kénnen, allgemein zu verstehen; sie kann unserer Auffas- 
sung nach héchstens in einzelnen Fallen fiir den konkreten Ablauf des Vereisungs- 
prozesses mitbestimmend und fordernd sein. Wir werden spater auf die von 
Montmory in Diskussion gezogenen Stoffe nochmals kurz zuriickkommen. 

Das Versagen einer umfassenden Erklarung des Vereisungsvorganges unter- 

kiihlter Wolken auf Grund kristalliner Vorstellungen legt es uns nahe, die expe- 

-rimentellen Befunde nach anderen Gesichtspunkten zu sichten, im Bestreben, 
allfallig zu einer allgemeineren Deutung der Eiskeimfahigkeit zu gelangen; doch 
muss vollstandigkeitshalber zunachst von zwei weiteren Erklarungsversuchen 
berichtet werden, die ebenfalls auf die Ergebnisse besonderer, teilweise sehr 
umfangreicher experimenteller Untersuchungen abstellen. 


B. Das Rausche Gefrierkernspektrum und die Struktur 
unterkuhlten Wassers 


In einer Grosszahl von Versuchen hat Rau [3] die Gefriertemperaturen von 
- Wassertropfen, die in einer geschlossenen Kaltekammer auf einer blanken 
Metallunterlage ruhen, bestimmt. Ordnet man seine Beobachtungen in der Weise, 


1) Physikalisches Institut der ETH. ; se . 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 493. 


~ 
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dass man die in einem kleinen Temperaturintervall ausgezahlten Erstarrunge: 
als Funktion der zugehorigen Gefriertemperatur darstellt, so erhalt man das nacl 
Rav benannte Gefrierkernspektrum, welches im Temperaturbereich von 0 bug 
— 30° drei Maxima zeigt, nadmlich bei — 4°, — 12° und — 19°. Wahrend ein oe 
stimmter Gefrierkern eines Tropfens bei Wiederholung des Versuches bisweilen 
eine Verminderung seiner Gefrierkeimfahigkeit zeigt, die sich in einem Absinken 
der Erstarrungstemperatur des Tropfens dussert, iibt die Alterung der Kerne au: 

die Lage der Maxima im Gefrierkernspektrum nur einen untergeordneten Einflus 

aus. So kénnen Versuchsbedingungen und Vorgeschichte der Kerne (Wasser- 
proben verschiedener Herkunft, direkt aus atmospharischem Niederschlag oder 
verschiedener Vorbehandlung usw.) sich in Gestalt und Héhe dieser Maxima in 
auffallender Weise, gelegentlich sogar bis zum vollstandigen Ausléschen eines 
Extremums, bemerkbar machen (wobei der Natur der Kerne nur eine sekundare 
Rolle zukommt), die Lage der Maxima bleibt aber in der Gesamttemperaturskala 
praktisch stets erhalten. Allerdings ist nicht klar zu erkennen, ob die Beobach- 
tungen von Rau, auf denen er sein Gefrierkernspektrum aufbaut, auch Messungen 
an eigentlichen Impfstoffen, wie AgI, mit einschliessen. Wenigstens sind tiber die 
Keimfahigkeit dieser Stoffe keine naheren Angaben zu finden. ; ; 

In der weiteren Entwicklung hat Rav versucht, einen plausiblen Zusammen- 
hang zwischen seinen Wahrnehmungen iiber die Unterkiihlbarkeit der Tropfen 
und den Laboratoriumsbeobachtungen hinsichtlich der Vereisung unterkiihlter 
Wolken herzustellen. Die von verschiedenen Meteorologen gefundene Haufig- 
keitskurve fiir die Temperatur von Schneewolken [4], [5], [6], [7] sind fiir Rau 
ein weiterer Beweis fiir die von ihm postulierte allgemeine Giiltigkeit seines 
Gefrierkernspektrums. : 

Wie wir im ersten Teil bereits erwahnt haben, zeigen die Einsatztemperaturen 
der von uns gemessenen eiskeimaktiven Stoffe Maxima bei Temperaturen von 
— 4 bis — 5° und von — 12 bis — 13°. Es ware aber nicht gerechtfertigt, diese 
Beobachtungen als weitere Stiitze, ja sogar als Erklarung fiir das Rausche Ge- 
frierkernspektrum heranziehen zu wollen; denn die Auslese der untersuchten 
Stoffe ist nach vorgezeichneten Richtlinien erfolgt, weshalb die Ergebnisse sta- 
tistisch nur bedingt im Spektrum eingeordnet werden diirfen. 

Wesentlich fiir unsere Uberlegungen ist, dass die von RAv verfochtene Allge- 
meingiiltigkeit des Gefrierkernspektrums ihn veranlasst, die Ursache fiir das 
Bestehen der 3 Maximalagen nicht mehr in Besonderheiten der Gefrierkerne, 
sondern einzig in strukturellen Verdnderungen des unterkiihlten Wassers zu 
suchen. Uber die Struktur des Wassers ist man im Temperaturgebiet tiber 0° seit 
den réntgenspektroskopischen Untersuchungen von BERNAL und FowLER [8] und 
weiterer neuerer Arbeiten hinreichend im Bild. So ist es gelungen, die Anomalie 
im Temperaturgang der Dichte des Wassers einwandfrei zu erklaren und auch die 
Druckabhangigkeit der Lage des Dichtemaximums zu verstehen. Nach den 
Beobachtungen von BERNAL und Fow er besitzt Wasser iiber 4°C eine vorwie- 
gend quarzartige Struktur, wahrend unter 4°C eine vorwiegend f-tridymitartige 
Struktur; in beiden Fallen ist der Aufbau durch eine tetraedrische Anordnung 
der O-Atome gekennzeichnet. Der Ubergang von einer Strukturart in die andere 
volizieht sich dabei allmahlich und stetig. Die Folgerungen von BERNAL und 
Fow Ler sind, abgesehen von der Lage der H-Atome, durch die spateren Beobach- 
tungen, insbesondere durch die rontgenanalytischen Arbeiten von MorGcAN und 
WaRREN [9] im grossen und ganzen bestatigt worden; allerdings musste die 
tridymitartige Struktur dem Begriff einer « zusammengebrochenen Eisstruktur » 
weichen. Es ist noch zu vermerken, dass gegensatzlich zu den Festk6rpern die 
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_Aussagen iiber die geometrische Anordnung der Atome in Fliissigkeiten nur auf 
Bre Entfernung zu Recht bestehen und vornehmlich statistische Bedeutung 
haben. 

Inzwischen sind die réntgenspektrographischen Messungen iiber die Struktur 
des fliissigen Wassers auch auf das Gebiet des unterkiihlten Zustandes ausgedehnt 
worden. Dorscu und Boyp [10] arbeiteten im Temperaturgebiet von + 21°C 
bis — 16°C, wobei ihre Réntgenlichtstreukurven auch im Gebiet unter 0° dieselbe 
allmahliche und stetige Strukturanderung des Wassers anzeigen, wie sie schon 
bei Temperaturen iiber dem Gefrierpunkt zur Geniige bekannt sind1), Ihre Auf- 
nahmen kénnen nur so verstanden werden, dass, je tiefer die Temperatur unter 
0° sinkt, die Wassermolekiile um so mehr dazu neigen, sich in einer weniger eng- 

_gepackten Struktur anzuordnen, bei welchem Vorgang die Struktur des Wassers 

’ desto eisahnlicher wird. Irgendwelche neue strukturelle Besonderheiten des unter- 
kiihlten Wassers treten dabei nicht in Erscheinung; auch der Durchgang durch 
den Gefrierpunkt vollzieht sich unbemerkbar. Mit diesen Ergebnissen stehen, 
nebenbei vermerkt, die Dichtemessungen von MoHLER [11] an unterkiihltem 
Wasser vollstandig im Einklang. Bei dieser Sachlage ist es nicht ohne weiteres 
verstandlich, wie weit die drei von Rau beobachteten Maxima des Gefrierkern- 
spektrums mit den oben erwahnten réntgenspektroskopischen Befunden ver- 
traglich sind. Eher diirfte infolge der kontinuierlichen Strukturanderungen ein 
einziges breites Maximum des Gefrierkernspektrums zu erwarten sein. Jedenfalls 
sind noch eingehende experimentelle Untersuchungen, die einen direkten Zu- 
sammenhang zwischen der Lage der Gefrierkernmaxima und den mutmasslichen 
Strukturbesonderheiten des Wassers aufzudecken hatten, vonnoten, wenn letztere 
fiir das Auftreten der Extrema verantwortlich gemacht werden sollen. Solange 
dies nicht der Fall ist, bleibt die Deutungsweise von Rau mehr oder weniger ein 
Postulat, trotzdem unverkennbar Griinde fiir die Heranziehung des Wassers zur 

Erklarung des Gefrierkernspektrumverlaufs bestehen. 


C. Adsorption von Wasserdampf an eiskeimaktiven Impfstoffen 
nach Birstein 


Kiirzlich hat BrrsTEIN [12] sehr sorgfaltig durchgefiihrte Messungen iiber die 
Adsorption von Wasserdampf an den beiden eiskeimaktiven Stoffen Silberjodid 
(AgI) und Bleijodid (PbI,) veréffentlicht. Diese im Temperaturbereich von 
4+ 26°C und — 20°C ausgefiihrten Untersuchungen haben gezeigt, dass die Ober- 
flache dieser beiden Impfstoffe solche Mengen an Wasserdampf aufnehmen kann, 
wie sie an gewohnlichen, nicht eiskeimfahigen Stoffen nie iiblich sind und daher 
als kennzeichnend fiir eiskeimaktive Stoffe vermutet werden miissen. Ferner 
zeigten bei einer Temperatur von — 20° unternommene Messungen einwandfrei, 
dass eine Adsorption von Wasserdampf an der Oberflache von Agl und PblI, 
nicht notwendigerweise Wassersdttigung zur Bedingung hat. Die beobachteten 
Adsorptionskurven nahern sich eindeutig und sanft den Eissdttigungswerten, 
wobei die Eisbildung auf der Oberflache der Impfstoffe bereits vor Erreichung 
dieser Sadttigung beginnt. Nach BrrstEIn ist demzufolge anzunehmen, dass 
bereits schon die ersten adsorbierten H,O-Molekiillagen eine eisahnliche Struktur 


1) Die markanteste Spitze der Streukurve (Intensitat als Funktion des Streuwinkels) weist eine 
stetige Verschiebung des Streuwinkels von 13°15’ bei 21°C bis 12°40’ bei=-16°C auf, wahrend das 
zweite Streumaximum. praktisch in unveranderter Lage bleibt. Die auffallendste Erscheinung ist 
dabei, wie sich das Minimum zwischen den beiden Hauptmaxima im unterkiihlten Zustand des 
Wassers rasch absenkt und das zweite Streumaximum an Scharfe entsprechend gewinnt. 
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besitzen. Weist in der Folge die entstandene Schicht eine geniigende Dicke (ver. 
mutlich von der Gréssenordnung von 200 Molekiildurchmessern) auf, kommt ihr 
Struktur schon ausgesprochen einem wahren Eiskristall nahe. | 

Die Ergebnisse und Folgerungen BirsTEIns stehen teilweise im Gegensatz z 
den von verschiedener Seite, so zum Beispiel von WEICKMANN [13]1) und AUFM 
KamPE[14] vertretenen Auffassung, dass fiir eine merkliche Adsorption von Wasser~ 
dampf an der Oberflache von Stoffen auch unter 0°C praktisch Wassersattigung 
notwendig ist und dass eine Sublimation, das heisst ein direkter Ubergang von 
dampfformigem Zustand in den festen, iiberhaupt nicht méglich ist. Jedenfalls- 
ist es BIRSTEIN gelungen, in einer Kiihltruhe, die vollstandig trépfchenfrei war 
und an den Wanden nur Eis enthielt, nach Ablauf von ungefahr einer Minute 
nach dem Einstauben von Silberjodid Eiskristallchen zu beobachten / 

Im Sinne der Erklarungsweise von BrrsTEIN stellt die Eiskeimfahigkeit nicht: 
léslicher Impfstoffe eine Oberflachenreaktion dar, die sich ausserordentlich emp- 
findlich auf den Zustand der Oberflache der festen Impfpartikel erweist. Er 
erachtet deshalb die Erweiterung seiner Adsorptionsmessungen auf durch Fremd- 
einwirkung gestorte Oberflachen als notwendig, um so weiteren Einblick in die 
Natur des Keimfahigkeitsmechanismus zu erhalten. Doch darf nicht iibersehen 
werden, dass BirsTEIN selber darauf hinweist, dass die Ergebnisse seiner Unter- 
suchungen nicht ausschliessen, dass gleichzeitig noch andere MOglichkeiten fiir 
eine Erklarung der Eiskeimbildung bestehen kénnen. 


D. Erweiterung der Weylschen Theorie der Gefrierkeimfahigkeit 


In. seiner Arbeit Surface Structure of Water and Some of Its Physical and 
Chemical Manifestations erbringt WryL [15] den Nachweis, dass die Oberflache 
fliissigen Wassers einen wesentlich anderen strukturellen Aufbau als das Innere 
besitzen muss. Die Wasserstoffatome, deren zugehérige Elektronen zufolge des 
Zustandekommens der chemischen Bindung praktisch in die Elektronenhiille des 
O-Atoms tibergetreten sind, weisen gegeniiber den Sauerstoffatomen eine nur noch 
verschwindende Polarisierbarkeit auf. Unter diesen Umstanden kann aber in der 
molekularen Anordnung der dussersten Oberflachenschichten Gleichgewicht nur 
dann bestehen, wenn die Wasserstoffkationen durchschnittlich eine tiefere Lage 
einnehmen als die Sauerstoffanionen, wodurch sich an der Oberflache eine elek- 
trische Doppelschicht ausbildet, deren negative Seite aussen liegt. Diese Erschei- 
nung kann auch so gedeutet werden, dass die den einzelnen Wassermolekiilen 
zukommenden elektrischen Dipolmomente fiir jene Molekiile, die nahe der Ober- 
flache sitzen, im Mittel nach innen gerichtet sind, Ein Wassertropfchen besteht 
demzufolge grosse modo aus einem innern Wasserkern mit quasi tetraedrischer 
Struktur, wie sie im Abschnitt B erértert wurde, und einer Oberflachenhaut, 
deren molekulare Bauelemente mit elektrischen Dipolmomenten ausgeriistet 
sind, die vorwiegend nach innen zeigen. 

Nach der Auffassung WeyLs liegt gerade im Bestehen einer Ausrichtung der 
molekularen Dipolmomente in der Oberflachenschicht eine der Hauptursachen fiir 
die leichte Unterkiihlbarkeit der Wassertropfen. Jedenfalls stellt die gleichsinnige 
Orientierung der Dipolmomente ein Hindernis fiir die Einleitung des Gefrier- 
prozesses von aussen (das heisst durch Warmeentzug) dar, indem sie den Uber- 
gang zur tetraedrischen Grundstruktur des Eises hemmt. Erst wenn durch irgend- 


1) Die Notiz WreickMANNs beschaftigt sich in erster Linie mit der Frage der Begiinstigung 
der Gefrierkeimfahigkeit durch kristalline Verwandtschaft 
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welche Einfliisse eine Abschwachung der Doppelschicht erfolgt und eine mehr 
oder weniger ausgepragte Entkopplung der ausgerichteten molekularen Dipole 
zustande gekommen ist, kann der Vereisungsvorgang des Tropfens bei einer ver- 
haltnismassig noch unbedeutenden Unterkiihlung seinen Anfang nehmen. 

In der Einfiihrung weist Wevt darauf hin, dass auch Ionenkristalle, wie zum 
Beispiel Natriumchlorid, bei vollstandig trockener Oberflache eine ahnliche 
Struktur mit zugehOriger elektrischer Doppelschicht besitzen miissen, weil auch 
hier dem Kation (Alkalimetall) eine merklich kleinere Polarisierbarkeit als dem 
Anion (Halogen) zukommt. Bei Zutritt von Feuchtigkeit wird aber, so folgern 
wir, unmittelbar der Lésungsvorgang eingeleitet und die besondere Oberflachen- 
struktur mit ihrer charakteristischen Doppelschicht geléscht. Auch bei unléslichen, 
vorwiegend homéopolaren Kristallverbindungen ist anzunehmen, dass unter 
“Umstanden die Oberflachenschicht in ahnlicher Weise wie beim NaCl mit einer 
elektrischen Doppelschicht ausgeriistet ist. Dies wird bei zweiatomiger Bindung 
zum Beispiel dann zu erwarten sein, wenn dem als Kation anzusprechenden Atom 
eine kleinere Polarisierbarkeit als seinem Partner zukommt. Auf dieses Bild 
stiitzend, hatte WEYL versucht, die Eiskeimfahigkeit von AgI durch die Annahme 
zu erklaren, dass ein AgI-Impfteilchen ebenfalls eine oberflachliche elektrische 
Doppelschicht (negative Seite aussen, wie bei der Wasseroberflache) besitzt, die 
beim Eindringen des Teilchens in die Wasseroberflache eine Auflockerung der 
Ausrichtung der molekularen Wasserdipole bewirkt. Nach unserer Auffassung 
ist es aber bei diesem Vorgang nicht unbedingt notwendig, dass das Impfteilchen 
tatsdchlich in die Oberflache des Wassertrépfchens eindringt, sondern es geniigt 
vollstandig, wenn dieses sich dem Wassertropfchen auf molekulare Entfernung 
nahert; denn die bei dieser Annaherung in Erscheinung tretende Wechselwirkung 
der Doppelschichten von Impfteilchen und Wassertropfchen, deren Momente im 
Annaherungspunkt in entgegengesetztem Sinn zeigen, kann bereits die Auflocke- 
rung der molekularen Dipolorientierung in der Oberflache des Wassertropfens 

‘nach sich ziehen. Méglicherweise ist damit die stets im Laboratorium beobachtete 
Ausserst rasche Vereisung unterkiihlter Wolken beim Impfen mit AgI noch ver- 
standlicher gemacht. 

Da dem Impfteilchen bei der eben gegebenen Erklarung des Vereisungsvor- 
ganges unterkiihlter Wolken im strengen Sinne die Rolle eines Gefrierkernes zu- 
kommt, sind wir damit zu einer molekularen Deutung der Gefrierkeimfahigkeit 
gelangt, die aber noch der experimentellen Untermauerung bedarf. Jedenfalls 
war bis dahin eine nur bescheidene Zahl von aktiven Gefrierkeimstoffen bekannt, 
die eine systematische Erweiterung der Weilschen Vorstellung noch nicht zuliess. 
Weil in der obigen Deutung die Gefrierkeimfahigkeit der Stoffe als reine 

' Oberflachenwirkung aufzufassen ist, schien es angezeigt, die von uns gefundenen 
Beobachtungsergebnisse tiber die Eiskeimfahigkeit der Stoffe im Rahmen des 
periodischen Systems der Elemente darzustellen. Figur 1 zeigt eine Tabelle des 
periodischen Systems, in welcher die Elemente gefrierkeimaktiver Stoffe mit- . 
einander durch eine Linie verbunden sind; dabei beschrankt sich die Wiedergabe 
auf Stoffe, deren Zusammensetzung nur zwei verschiedenartige Elemente enthalt. 
Stoffe, die schon im Temperaturgebiet von — 10° bis 0° sich als aktiv erweisen, 
sind durch eine kraftig ausgezogene Linie kenntlich gemacht. Alle Linien beginnen 
grundsatzlich bei einem Element auf der linken Seite des periodischen Systems 
und fiibren mit den entsprechenden Verzweigungen nach ihren Partnern, die 
stets auf der rechten Seite des Systems liegen. In der Tabeile sind zudem die 
Werte der Polarisietbarkeiten (« - 1025 cm—3), welche die einzelnen Elemente als 
freie Ionen besitzen, angegeben ; die mit einem Kreuz gekennzeichneten Zahlen 


Figur 1 


Gefrierkeimfahigkett der Stoffe und Periodisches System 
™= aktiv zwischen 0°C und — 10°C. 
— aktiv zwischen — 10°C und — 20°C. 
Zahlenangaben: Polarisierbarkeiten der freien Ionen (« +1025 cm 3). 


sind vermittels den von L. PAULING angegebenen Werten fiir die lonenrefraktionen 
berechnet. 

Wir entnehmen der Tabelle auf den ersten Blick, dass jene Stoffe, die als 
eiskeimaktiv befunden worden sind, ohne Ausnahme sich stets aus solchen Atomen 
zusammensetzen, deren Polarisierbarkeiten sich merklich unterscheiden. Dabei iiber- 
nimmt stets das Element mit der kleineren Polarisierbarkeit entsprechend seiner 
linksseitigen Lage im peviodischen System mehr oder weniger betont die Rolle eines 
Kations, wahvend seinem Partner jene des Anions zukommt. Im Umstand aber, 
dass alle gefrierkeimfahig befundenen Stoffe in Wasser nur schwer léslich sind}), 
erkennen wir einen Hinweis dafiir, dass die chemische Bindung schon merklich 
homéopolaren Charakters ist. Klar und deutlich bringt die Tabelle zum Ausdruck, 
dass fiir die Gefrierkeimfahigkeit nicht die kristallinen Eigenschaften eines Impf- 
stoffes verantwortlich sind, sondern vorwiegend jene molekularen Krafte, die die 
Struktuyr der Oberfldche der Stoffe beherrschen. 

Da voraussichtlich die elektrische Doppelschicht der eiskeimaktiven Stoffe 
um so ausgepragter sein wird, je starker die Polarisierbarkeiten von Kation und 
Anion voneinander abweichen, werden wir vermuten, dass bei festgehaltenem 
Kation und variiertem Anion die Gefrierkeimfahigkeit der Stoffe mit wachsender 
Polarisierbarkeit des Anions zunimmt. Tatsachlich haben wir im ersten Teil 
dieser Arbeit feststellen kénnen, dass unsere Vermutung im grossen und ganzen 


1) Auf den Umstand, dass gefrierkeimaktive Stoffe wasserunldslich sind, hat schon H. K. 
WEICKMANN [13], S. 313, hingewiesen. 


Vol. VI, 1955 Kurze Mitteilungen - Brief Reports - Communications bréves 491 


stimmt, unbeschadet, ob es sich dabei um zwei- oder mehratomige Impfstoffe 
handelt. Dieses Vorgehen ist aber nicht streng korrekt; denn die durch die 
chemische Bindung vorgezeichnete Elektronenkonfiguration muss fiir die Impf- 
stoffe, deren Gefrierkeimfahigkeit gegeneinander abgewogen wird, einen 
gleichartigen Aufbau besitzen. Wiirden wir zum Beispiel ungeachtet der chemi- 
schen Bindung bei konstant gehaltenem Anion und variiertem Kation den nam- 
lichen Vergleich wie oben im umgekehrten Fall ziehen, so kénnten wir tiberhaupt 
keinen unmittelbaren Zusammenhang zwischen Gefrierkeimfahigkeit und Polari- 
sierbarkeit aus unseren Beobachtungen erkennen. 

Vergleichsreihen, von denen wir unter Umstanden erwarten diirfen, dass sie 
einen Einblick in das Verhaltnis zwischen Polarisierbarkeit und Keimfahigkeit 
von zwei- und allfallig auch mehratomigen Impfstoffen vermitteln, erhalten wir, 
wenn wir zum Beispiel bei festgehaltenem kationartigem bzw. anionartigem 
Element die Partnerelemente ein und derselben rechts bzw. linksseitigen Kolonne 
(Unterkolonne) des periodischen Systems entnehmen. Als Vertreter dafiir k6nnen 
wir die folgenden Reihen versuchsweise herausgreifen: 


a) mit festgehaltenem kationartigem Element 


ie CuCl CuBr Cul?) 
n n az 
2; AgCl AgBr Agl 
£17 f —4,0 
3 CuO CuS (Cusis CuTe 
mete? —4,0 a4 } 
4, ZnO ZnS ZnSe Cole 
f n 12,3 if 
De Ag,O Ages Agese AgoTe 
; —6,2 mae its f 
6. HgCl, HgBr, Hel, 
n —13,0 
b) mit festgehalienem anionartigem Element 
ths CuCl AgCl AuCl 
n —11,0 ij 
8 Cul Agl Aul 
—4,2 —5,0 —8,0 
9, ZNSE CdSe HgSe 
S123 —4,2 f 
10. Cu,O Ag,O Au,O 
—5,1 —6,2 i 
dal; Hel, Cdl, Znl, 
1x0 +15,0 f 


Die aufgefiihrten Ver 
unsere Erwartung iiber die Abhangigkeit 
tisierbarkeit des Anions im Falle der zw 
Falle der dreiatomigen nur in der Versuc 


1) Die angefiihrten Zahlen bedeu 
Buchstaben 7: nicht aktiv bis — 18°; un 


suchsreihen zeigen, dass bei festgehaltenem Kation 
der Gefrierkeimfahigkeit von der Pola- 
eiatomigen Verbindungen zutrifft, im 
hsreihe 6, wahrend °5 sich gegenlaufig 


ten die Einsatztemperatur der Gefrierkeimfahigkeit; die 
d f: Beobachtungen fehlen. 


~~ 
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erweist. Bei festgehaltenem Anion hingegen zeigt sich im Falle der zweiatomigi 
Verbindungen nur die Versuchsreihe 8 unseren Erwartungen gemass (Abnahme 
der Keimfahigkeit mit zunehmender Polarisierbarkeit des Kations). Die beiden 
Reihen 7 und 9 verhalten sich gegenlaufig. Die fiir den Fall von dreiatomiger 
Verbindung aufgezahlten Versuchsreihen 10 und 11 zeigen sich, soweit wir die 
Sache bei den sparlichen Angaben iiberblicken koénnen, praktisch unberiihrt von 
der Anderung der Polarisierbarkeit des Kations. ‘ 

Das Beobachtungsmaterial weist noch zu viele Liicken auf, um ein abschlies- 
sendes Urteil tiber die hier aufgeworfene Frage des Einflusses der Veranderlich- 
keit der Polarisierbarkeit auf die Gefrierkeimfahigkeit zu fallen. Von den fiir eine 
Beurteilung wichtigen zweiatomigen Impfstoffen fehlen unter anderem das leicht 
zersetzbare Silberbromid (AgBr) und die Telluride von Cu, Zn und Ag. Auch 
ware es bei der Sparlichkeit an Messungen noch verfriiht, weitere Versuchsreihen 
aufstellen und erértern zu wollen. Vielmehr ist es unumganglich, das Erfahrungs- 
material vorerst merklich zu erganzen. Schliesslich ist an dieser Stelle noch auf 
die von Montmory hinsichtlich des Epitaxieeinflusses in Diskussion gezogenen 
Manganverbindungen Mg,Sb, und MgTe zuriickzukommen, um festzustellen, dass} 
diese beiden als eiskeimaktiv erkannten Stoffe sich ohne weiteres in das von uns) 
in Tabelle 1 aufgestellte Schema einfiigen lassen}). | 

Fiir die Klarung der Frage nach der Abhangigkeit der Eiskeimfahigkeit der’ 
Stoffe von den Polarisierbarkeiten der beteiligten Atome, oder allgemein aus-. 
gedriickt, von den Oberflachenkraften zwischen Impfteilchen und Wassertrépf-- 
chen, werden Versuche iiber die Gefrierkeimfahigkeit, bei denen die Oberflache: 
der Impfteilchen durch kiinstliche Beeinflussung gestért wird, besonders wertvoll, 
wenn nicht sogar entscheidend sein. Vorteilhafterweise hatten sich diese Versuche ' 
zunachst auf zweiatomige Stoffe zu beschranken. Liegen die Ergebnisse solcher’ 
Messungen vor und sind die vorhin erwahnten Liicken im Erfahrungsmaterial. 
grosstenteils geschlossen, wird auch der Zeitpunkt gekommen sein, sich mit dem. 
Beobachtungsmaterial von Rav und insbesondere mit der von ihm vertretenen. 
Allgemeingiiltigkeit des Gefrierkeimspektrums auseinanderzusetzen. | 

Es darf nicht tibersehen werden, dass die bei gefrierkeimaktiven Stoffen ge-. 
fundene starke Unterschiedlichkeit der Polarisierbarkeiten nicht als eine hin-. 
reichende Bedingung fiir das Auftreten der Gefrierkeimfahigkeit betrachtet wer-. 
den darf. Vermutlich werden Stoffe bestehen, die dieser Bedingung gentigen und 
dennoch nicht gefrierkeimaktiv sind. Eine zweite Bedingung scheint erfahrungs- | 
gemass ebenfalls erfiillt sein zu miissen, ndmlich die Unléslichkeit (oder zumin- 
destens nur dusserst geringe Léslichkeit) des Impfstoffes in Wasser, wodurch die 
Léschung der oberflachlichen Doppelschicht durch den Loésungsvorgang nicht in 
Erscheinung treten kann. Dies deutet darauf hin, dass lésliche Impfstoffe nur 
bei vollstandig trockener Oberflache im Sinne der erweiterten Weylschen Vorstel- 
lungen gefrierkeimaktiv sein kénnten. Dem steht aber entgegen, dass ldsliche 
Stoffe allgemein eine merkliche Hygroskopizitat zeigen. 

Abschliessend muss noch ausgesagt werden, dass dem in diesem Abschnitt 
erorterten Versuch einer Erklarung der Gefrierkeimfahigkeit schlechthin nur. 
die Bedeutung eines Diskussionsvorschlages zukommt. Ihm mehr beizumessen, 
ware in Anbetracht der, gesamthaft betrachtet, noch merklichen Sparlichkeit des 
Beobachtungsmaterials zu gewagt. Auch ist der Gedanke BIRSTEINS, dass unter 

1) S. S. Brrsrein und C. E. ANDERSEN [16] ver6ffentlichten kiirzlich Messungen tiber die Ge- | 
frierkeimfahigkeit von Stoffen unter gleichzeitiger Bezugnahme auf friihere Messungen von C. L. 
Hosrer [17]. Da alle Beobachtungen sich auf die Temperatur von — 20°C beschranken, nehmen | 
wir von einer direkten Gegentiberstellung mit den von uns erérterten Messergebnissen Umgang. 
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Umstanden nicht nur ein einziger Mechanismus fiir die Eiskeimfahigkeit verant- 
wortlich gemacht werden kann, ebenfalls noch weiter zu verfolgen. 
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Summary 


The results of the measurements published in the first part of this report are 
compared to the existant attempts at explanation of the ice nucleability of 
substances. Of these only the interpretation of the nucleability of supercooled 

clouds proposed by WEYL can, somewhat extended, be deemed to fit actual 
observations. It is found in the case of 2-atomic seeding substances that one of 
the elements (the one on the left-hand side of the periodic system) must be 
remarkably less polarizable than the partner (on the right-hand side of the 

_ periodic system). Another characteristic property of the seeding substances found 
to be active is their non-solubility in water. The forces responsible for the 
nucleability of the substances are those which also determine the structure of 
the surface of the seeding particles and of the water droplets. 


(Eingegangen: 6, September 1955.) 
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Zur Berechnung der Geschossabweichung unter dem Einfluss eines 
Seitenwindes 


Von Ernst ROTH-DESMEULES, Luzern 


* 
Erfolgt die Bewegung eines artilleristischen Geschosses nicht in einer ruhenden ¢ 
Atmosphare, sondern in einer senkrecht zur ungestorten Flugbahnebene bewegten 
(herrscht also Seitenwind), so fiihrt dies zu einer seitlichen Abweichung des Ge- 
schosses; diese kann nach der Naherungsformel : 
t 
} 


— 
= 
~— 


x 
he ee (7 U9 COS Hy 


berechnet werden. Hier bedeuten: z,, seitliche Abweichung, w, Windstarke, T | 
Flugzeit, X horizontale Flugdistanz (der ungestérten Flugbahn), v, Anfangsge- | 
schwindigkeit, #, Abschusswinkel. In den meisten Hand- und Lehrbiichern der | 
Ballistik wird hemenlce dass diese Formel auf Grund der vorgenommenen Ver- — 
nachlassigungen nur fiir steile Flugbahnen nicht mehr geniigend genau sei!). Nun 
zeigt aber die numerische Rechnung, dass im allgemeinen die Naherung auch fiir | 
sehr steile Bahnen, die nur wenige Grad vom Schuss senkrecht aufwarts ab- | 
weichen, noch ausserordentlich gut ist, indem der relative Fehler unter 10-4 
bleibt. Dies soll im folgenden naher untersucht werden. 

Es werde ein Koordinatensystem eingefiihrt, so dass die *-Achse horizontal 
in die Schussrichtung zeigt, die y-Achse senkrecht nach oben, und schliesslich — 
soll die z-Achse das rechtshandige System vervollstandigen. Betrachten wir zu-~ 
nachst die Flugbahn bei ruhender Luft, dann lauten die Differentialgleichungen — 
der in einer Ebene erfolgenden Geschossbewegung 


Xn =F =, in , ae re = @,, Vn 7 &(Yn) : (2) 


Anfangsbedingungen fiir ¢= 0: %, = Vn = 0, 4, = Ug COS Oo, Yn = Vy Sin B. Die 
Funktion ©, = ©, (4), Va» Vn) ist “proportional ‘dena durch die Gesch windiest 
Vy_ = (42 + ai? atte ‘dividierten Luftwiderstand, Der Index n deutet darauf hin, 
dass es sich um die Noyvmalbahn handelt. 
Weht nun ein Wind w, norinal zur Flugbahnebene (x, vy), dann nehmen die 
Bewegungsdifferentialgleichungen die Form an: 
ty = — ®,, Fry ’ Vw ae —D, Vw a § (Vw) » fy = = 2, (2 ro Wz) ' (3) 


Anfangsbedingungen fiir ¢ = 0: 


%y(0) Vw(0) = 2,(0) = 0, % w(0) = Up COS By = Xo, Vw(0) = V9 Sin Io = Vor Zy(0) = 
Die Geschwindigkeit ist v2 = 42 + y2 4 (2, — w,)%. Dies ist die sogenannte 
Windbahn. 

Vermoge der Transformation 


xy = %,COSY, Vy = Vp, Zy = Wt, — x, siny (4) 
1) Zum Beispiel: C.CRANz, Aussere Ballistik, Bd. 1 (Springer, Berlin 1925), S. 301. H. ArHEN | 
Ballistik (Quelle und Meyer, Leipzig 1941), S. 215. 


i! 
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werde ein mit der Atmosphare bewegtes Koordinatensystem eingefiihrt. Es ist 
ke cos afl 
mec a oe 


(5) 


Mf =a : v. 
. tgy ry Se: siny = (42 = ware : 
 Setzt man die Gréssen (4) in (3) ein, dann fallen die erste und dritte Differential- 
_ gleichung zusammen, und fiir die Relativbahn, die wiederum eine ebene Flugbahn 
ist, gilt 


K, = —O, fn Ve= —O, Ve — BV r) - ‘ (6) 
- Anfangsbedingungen fiir ¢= 0: 
; 4 (0 s 3 : 
, x= I= 0, 4,(0) = 2D = (28+ WHI, _9,(0) = Fo. 


In diesem Koordinatensystem haben wir also eine normale Flugbahn vor uns, 
allerdings mit neuen Anfangsbedingungen. Fiir die Anderung der Anfangselemente 
gilt in erster Naherung 

eo eee ae Cw 
sual ag FTE Ath= —Z t) tg dy. 
Die Projektion dieser Relativbahn auf die (¥,, ¥,)-Ebene gentgt denselben An- 


- fangsbedingungen wie die Normalbahn. Die Differentialgleichungen (2) und (6) 
unterscheiden sich jedoch im Werte der Funktion ©. Es ist 


‘ 


Dy es Dy (Uy, Yn) mit U;, = 4¢ + vo 


~ und 
re j —_ 2 2 
P= D,(v,, y,) Mit vi= UZ + wz - 
Ferner ist 
; W, 
rs Be = 0, = 0, COS, MIC Re rat 


0 


Wenn nun fiir t= 0 auch @, = ®, ist, dann stimmen nach dem Eindeutig- 
keitssatz fiir die Lésungen von Differentialgleichungen die Integrale im ganzen 
Verlauf iiberein. Nun ist, wenn nur die wesentlichen Grdssen hervorgehoben 


werden 
h 
ov, ») =“ Fw. 
_ Die obige Gleichheit ist damit erfiillt, wenn 
il 1 
as i (D),) = a3 F(v,) 
oder 
F(v, cos w) = F(v,) cos@ . 
Diese Beziehung ist nur méglich, wenn F eine lineare Funktion von v ist. 


F(v) =cv. (7) 


In diesem Falle wird die. Bewegung in der (¥, y)-Ebene unabhangig von der seit- 


lichen Abweichung . Ly fetes 
In Wirklichkeit ist die Funktion F(v) nicht allein von der Geschwindigkeit v, 


sondern von der Machschen Zahl M = v/s abhangig. Der Einfluss der von der 
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Hohe y (das heisst der Temperatur) abhangigen Schallgeschwindigkeit spielt 
hier aber keine grosse Rolle, und F ist, als Funktion von v allein betrachtet, 
nahezu linear, wie die Abbildung fiir zwei verschiedene Falle zeigt1). Dies gilt | 
im allgemeinen fiir Geschwindigkeiten, die grésser sind als die ephallgeecS an | 
keit, und zwar ebensogut fiir Flachbahnen wie auch fiir Steilbahnen. 4 
Fiir die seitliche Abweichung des Geschosses folgt aus den Transformations-— 


formeln (4) mit Hilfe von (5) : 
3 Ky W it 
Zy = wW,t— #, siny = w,t a , 
Xo : 
Da in unserem Falle unter der Voraussetzung (7) *,, = *, ist, so erhalten wir die~ 
in der Einleitung gegebene, nun exakt giiltige Formel (1) ‘ 
2 
sd i 
ay unl ph 


4 


wo auf der rechten Seite die Werte der Normalbahn (2) stehen. 

Bemerkenswert ist noch, dass sich die Beziehung +, ~ %, bereits aus der 
Kenntnis der Flugbahngréssen der Normalbahn — namlich aus der approxima- 
tiven Giiltigkeit der Beziehung (7) — ergibt, unabhangig von der zusatzlichen Be- 


| 
| 
| 

] 
} 


Figur 1 


Die Funktion /’(v) fiir zwei verschiedene Luftwiderstandsgesetze. 1 Flachbahn, 2 Steilbahn. 


stimmung der Einflusskoeffizienten fiir eine Avy- und eine 4% )-Stérung, die sonst 
erforderlich ware, um die Elemente der Relativbahn (6) aus der Normalbahn (2) 
zu berechnen. 

Die Formel (1) gibt uns den Einfluss eines konstanten, langs der ganzen 
Flugdauer wirkenden Seitenwindes. Es ist offensichtlich, wie sie umzuformen 
ist, wenn der Wind nicht zur Zeit ¢ = 0, sondern erst zu einem spatern Zeitpunkt 
einsetzt. Das Entsprechende gilt auch in diesem Falle. 


1) Eine additive Konstante ist, da cos @ nahe 1 liegt, praktisch belanglos. 
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-* Zusammenfassend kann damit gesagt werden, dass bei Uberschallgeschwin- 
_digkeiten die einfache Formel (1) den Einfluss des Seitenwindes sowohl fir 
flache als auch fiir steile Flugbahnen sehr gut wiedergibt. 


f 


Summary 


The first order correction which should be made to the z-coordinate of the 
standard trajectory of a projectile to account for a cross-wind is given by (1). It 
is shown that this formula is a very good approximation also for angles of depar- 
ture near 90°. 


(Eingegangen: 13. April 1955.) 


Bemerkungen zur numerischen Integration gewohnlicher 
Differentialgleichungen n-ter Ordnung 


Von HEINZ RUTISHAUSER, Ziirich?) 


Es ist bekannt, dass man eine Differentialgleichung »-ter Ordnung in ein 
System von 7 Differentialgleichungen erster Ordnung verwandeln kann; davon 
macht man insbesondere fiir die numerische Integration gerne Gebrauch. Es 
besteht aber auch die Moéglichkeit, eine Differentialgleichung »-ter Ordnung 
direkt — das heisst ohne Verwandlung in ein System — numerisch zu integrieren, 
namlich mit den Formeln von Kutta-NystrOmM?), bzw. Kutrra-ZURMUHL?). 

Abgesehen von anderen Vorteilen, wird den beiden letztgenannten Integra- 

-tionsverfahren eine auffallend grosse Genauigkeit nachgeriihmt. Der Verfasser 
-hat daher versucht, die verschiedenen Verfahren auf ihre Genauigkeit: zu ver- 
gleichen, und zwar durch Anwendung auf lineare Differentialgleichungen mit 
konstanten Koeffizienten. Fiir diese kann man namlich sowohl die Differential- 
gleichung als auch die durch das numerische Verfahren induzierte Differenzen- 

_ gleichung exakt lésen. Dadurch war es méglich, die Untersuchung zum grossten 
Teil analytisch durchzufiihren, einige zusatzliche numerische Rechnungen wurden 
mit der Z4 durchgefiihrt 

Die Herleitung der Ergebnisse ist in einem internen Bericht des Instituts fiir 
angewandte Mathematik der ETH. festgehalten; wir geben hier lediglich die 

- Resultate wieder. 

1. Die Formeln von RunceE-Kutra (fiir Differentialgleichungen 1. Ordnung), 
Kutra-Nystroo (fiir Differentialgleichungen 2. Ordnung) und Kutra-ZURMUHL 
(fiir Differentialgleichungen n-ter Ordnung, » > 2) haben alle dieselb Genauig- 
‘keitsordnung 4, das heisst, dev Fehler bei numerischer Integration mit der Schritt- 
weite h tiber ein festes Intervall ist bei allen diesen Verfahren ~c h4. Dies gilt auch 
dann, wenn die Differentialgleichung von der Form y) = f(%, y) ist; allerdings 
wird dann der Koeffizient c (vor h*) mit wachsendem » kleiner. Die verschiedent- 
lich behauptete Vergrosserung der Genauigkeitsordnung mit wachsendem » trifft 
nur fiir die ersten Integrationsschritte zu und ist darum ohne praktische Bedeu- 


1) Institut fiir angewandte Mathematik der ETH. 


2) Siehe Literatur [1] oder auch [4], Kap. Besides 
3) Siehe Literatur [2] oder auch [4], Kap. I, § 3.2. 


ZAMP VI1/32 
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tung. Uberdies hat der Verfasser festgestellt, dass bei Integration der Gleichun 


NW 6 
(4. z 1)" eee (") yD (") cn heme canbe VES, 


nach KutTa-ZuRMUHL auffallend grosse Fehler entstehen. 

2. In ahnlicher Weise untersuchte der Verfasser die Integrationsformeln von _ 
ApAms-FaLKNeER!), die man durch Ubertragung der Adamschen Methode auf 
Differentialgleichungen hdherer Ordnung erhalt. Die Resultate dieser Unter-_ 
suchungen sind den unter 1. genannten Ergebnissen véllig analog. Die Genauig- 
keitsovdnung des Verfahvens von ADAMS-FALKNER (Extrapolation oder Interpola- 
tion) hdéngt nur von der Anzahl der mitgefihrten Differenzen A®f in den Formeln 
(5. 22) bzw. (5. 32) bet CotLatz [4] Kap. I, ab, jedoch nicht von n. - 
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Summary 


The author reports on the results of an investigation concerning the accuracy — 
of the so called Kutta-Zurmiihl formula [2] for numerical integration of ordinary | 
n-th order differential equations. If one integrates with step-length h over a_ 
fixed interval, the error of this formula is ~ h* (for small h), irrespective of 1 
(and not ~ h"*+* as it has been claimed). ; 


(Eingegangen: 14. Mai 1955.) 


. nme 


Torsion of a Circular Shaft Press-fitted with a Disc 


By HrIbEo Sarto, Sendai, Japan?) 


For a shaft of varying circular cross-section under the action of torsion, it is | 
shown that the displacement of any point is directed at right angles to the axial | 
plane passing through the point’). If this displacement is denoted by v and if’ 
polar coordinates (r, 0, z) are used with the axis of the shaft as the Z-axis, the © 
stress distribution in the shaft is determined by a solution of the partial differen- — 
tial equation 

0? ioe ou @) 0% 
parctapss Og) pp (2) 


1) Siehe Literatur [3], oder auch [4], Kap. I, § 5.4, und 5.5. | 
2) Department of Mechanical Engineering. The Tohoku University. | 
3) S. TrmosHEenko, Theory of Elasticity (McGraw-Hill Book Comp., 1934), py 226. 
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and the stresses are given by 


ou 

¢ iy 5 m=Gr5-(2), (2) 

_ where G is the modulus of rigidity. In this report, we discuss the torsional problem 
for a circular shaft press-fitted with a disc or a wheel. The outer diameters of the 
shaft and the disc are denoted by 2a and 25, the thickness of the disc is # (Figure 1). 
In order to simplify the analysis, we shall proceed with our calculations assuming 
that a slip does not occur on the contact surface between the shaft and the disc 
and that at a great distance from the disc the stress in the shaft becomes a state 
of pure shear and its moment is M. For the shaft and the disc, the notations I 

and II are used respectively. 


Figure 1 


Circular shaft press-fitted with a disc 


The boundary conditions to be satisfied are 


(a) ci=0 (¢=0, h), (3) 
h 
230 b* t,dz=—M (7—= 0). (4) 
0 
(b) On 7 = a 
t= Ty (0 S2°< A), t= 0 (z>h), (5) 
eee oi, Sr any. (6) 
(c) i 
a a 4) (Z = oo) f (7) 


where t is the circumferential shearing stress of the simply twisted part. 
Considering the boundary conditions, we may introduce the following expres- 

sions for displacements and stresses. 

For the disc: 


2 cS Ky(%, ”) 
a oe a Bake ae 
Gu —T |4o ay +e Ay ACE COS &, | 
2 SS K,(a ; 
I aon 2\%n z\, 8 
co |4, 7? : x ai Ko(a%, & an | 


ee IRGC, 9) nes ( M at) 
re oe Mie EAS z A Ne 
Vo ie aya Ca te ie 
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For the shaft: 


cas | r+ 38 fry 2 ee | 608 a # das fs 


co 
Mi Ginitiss [Fe dies cosa * da, 


no -28 fry 


_ where J and K are modified Bessel functions of the first and the second kinds, 
A,, B are arbitrary constants, and F(«) is a function of «. 


tae already satisfies the condition (3). Also onr=a 


Il = 
T,g=T Agrees Ay, CO8E es & : 
n= 


- From condition (5) and a Fourier transformation follows 


sin (« h/a) 


+34 (-1)" a sin (a h/a) 


(Gt) eA 7 i eee ay AE 


Substituting this into as (9) and equating v! to v!!, we have by the usual 
Fourier analysis 1 


: 
: 
: 
: 
t 
; 


I 5 % ee 
Mandeep) allt pete, 2 dabei ae | 


where 


anf ) sin? ae =f xh a) sin? (« A/a) d 


a se se 
for) (11) 
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The integrals of equation (11) follow from the theorem of residue: 


hla oo 1 — g-4m2hla 


co 
Ing = — ZS na (~ + s), (12) 
1 TN Pa Te (Cer c) ie 
ik Galen ear Sapte ack = 8) 


where 
a (1 ial e 4m 2 hlay 
mS (AR +e, 2?) (AR + OX, 2?) 


Me 


(13) 


and 4/,, is the m-th root of 4-2 J,(4) = 0, in which J, is the Bessel function of the 
first kind. an 
Substituting equation (12) into equation (10), we have 


lee) 
A,— EAn(-1U"Qne=Aois (8= 1,2, 3,-..), (14) 
n= 
where 
Gl (ay? (—1)8s Ky(a,4). & 
= m 
Qns 13 G (;;) 6, Ky (a, @) J ’ns 
’ Gt pat? (a1) 62K .(a,.2) ; 
4 Ss Cae} 15 
Meron dhrecica em 5) 


A: G") K,(a, a) Ty(a, 4) 
Sima e fie eile 


Now, we shall determine the value of A). The limiting value of ee in equation 
(9) for ¥ = oo becomes 
lim th. = —4T 4(4) ee 


X00 


Comparing this with condition (7), we find 


and condition (4) also is satisfied. 
The series of d,,, do not converge rapidly when » and s are large. For large 


values of m, we therefore substitute for /,, its asymptotic value [m + (3/4) 7]). 
Applying EuLER’s summation formula?) to the series directly, we have 


lee) m’ --1 al 48 Ae (x a)? 
m et m ihe fesce ee 
ae, ems = mst Trak (ak ad) 8 28, + (ay a)? 
i" Ani LS 
2 (Ae, + Gp G2) (AB, + Oy a) 


If we take the appropriate value of m’, the sum of the series is precisely calculated. 


ay WATSON, Theory of Bessel Functions (Cambridge University Press, 1922), p. 507. 
<3 1. A. Bromwicn, Theory of Infinite Series (Macmillan & Co., London, 1926), p. 304. 


G.N. 
‘oe 


~~ 
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The unknown constants A, can be obtained in each particular case from 
equation (14) by the method of successive approximations. Numerical calculations 
have been carried out for the case h/a = 2/2, b/a = 2, G! = G"!. Taking the first 11 
terms of the series, we obtain the following values of 4,,: 


Table 
a ee renee eee 
An 0:2281 | — 0:1536 0:1147 | — 0-:0921 0-0774 | — 0:0672 
[Te | 7 | 8 | 9 | 10 | ll 
Zale 0:0595 | — 0-0536 00-0489 | — 0°0450 0-0417 | 


Actually, the convergency is very slow for stresses on y = a. Hence methods 
similar to those used in a previous paper!) have been applied for the computation 
of the stresses on y = a. In Figures 2 and 3, the stress distributions in the shaft are 


shown. The curves of constant resultant shearing stress (/x3, + 1§,/t = constant) 
and v/y = constant are shown in Figure 4. 


F 
se a at et ln patel piglet ait = 9h penn en wise nl eens aaapas 


Wes, 


0 ee ee ; 
D9 02 O04 O6 08 10 12 14 16 18 O 02 O04 08 O08 10 72 14h) 
Zh —»> 2/h ie | 

Figure 2 Figure 3 
Distribution of Ei in the shaft. Distribution of t in theshaft. | 


1) H. Sarro, Techn. Rep. Tohoku University 18, No. 2, 178 (1954). 


Vol. VI, 1955 Varia — Miscellaneous — Divers : 503 


Figure 4 
—— Curves of yz 2+72 /¢=constant; ————: Curves of v/r = constant. 
Zusammenfassung 


Es wird hier ein Lésungsverfahren besprochen fiir das Torsionsproblem der 
Welle, auf die eine Scheibe oder ein Rad gepresst ist. 


(Received: April 13, 1955.) 
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Herbsttagung der Schweizerischen Physikalischen Gesellschaft 
vom 24. bis 25. September 1955 in Pruntrut 


Untersetzer mit Flachentransistoren'), von E. BaLpINGER und M. NICOLET, 
Basel?). 

Untersetzer mit Flachentransistoren sind schon verschiedentlich beschrieben 
worden [1], [2], [3], [4], [5], [6], [9]*), doch waren Angaben iiber das Verhalten ein 
und derselben Schaltung unter verschiedenen Betriebsbedingungen nicht erhaltlich. 
Deshalb ist der hier beschriebene Untersetzer sowohl auf seine Temperaturab- 
hangigkeit als auch auf sein Auflésungsvermogen bel regelmassig und bei stati- 
stisch eintreffenden Eingangsimpulsen untersucht worden. 


1) Vorgetragen an der Herbsttagung der Schweizerischen Physikalischen Gesellschaft, mit 
Ausnahme der Ergebnisse, die sich auf Philco-Transistoren beziehen. 


2) Physikalisches Institut der Universitat Basel. 
3) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 506. 
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Das in Figur 1 dargestellte Schema des Untersetzers entspricht demjenigen 
eines iiblichen Triodenscalers mit zwei stabilen Zustanden. Bei Schaltungen 
dieser Art ist es wesentlich, den Arbeitspunkt des leitenden Transistors so zu 
wahlen, dass die Basis-Kollektorspannung einen gewissen Minimalwert nie unter- 
schreitet, ansonst das Impulsverhalten des Transistors beeintrachtigt und somit 


Figur 1 
Schaltschema des Untersetzers mit Transistoren. 

Transistorentyp Philips OC 70 und 71 Philco SB-100 ; 
R, 15 kQ 15 kQ . 
re 3 kQ 3 kQ 
Rp 15 kQ hse vase ] 
Rs 3 RQ 3 kQ 
Rc 2 RQ 2 RQ 
Cy 1000 pF 20 pF | 
c. 400 pF 10 pF 
Cp 2,5 uF 2,5 uF 
Cr 2,5 uF 2,5 uF | 
L 2 mH 40 wH Sy 


das Auflésungsvermégen verschlechtert wird [7], [8]. Der Untersetzer spricht auf 
positive Eingangssignale an, welche den leitenden Transistor sperren. Das Aus- 
gangssignal wird zuerst differenziert, worauf eine Diode die negativen Impulse 
unterdriickt. Ist ein geringeres Auflésungsvermogen zulassig, so kénnen die 
einzelnen Stufen direkt mit kleinen Seriewiderstanden gekoppelt werden}). 

Aus allgemeinen Betrachtungen iiber den Ein- und Ausschaltvorgang bei 
Flachentransistoren [8] lasst sich folgern, dass das zeitliche Auflésungsvermégen 
von Kippschaltungen eine obere Schranke nicht iiberschreiten kann. Eine 
Schatzung zeigt, dass diese Schranke fiir Untersetzer unserer Art bei rund fo/2 
hegt (f) = Grenzfrequenz des Transistors). Das wirklich erreichbare Auflésungs- 


1) Auch eine induktive Kopplung tber die Kollektorinduktivitaten ware moglich, ist aber 
nicht untersucht worden. 
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vermogen hangt von der Schaltung ab, wobei spezielle Massnahmen — zum 
Beispiel das Einfitthren von Induktivitaéten im Kollektorkreis (L, Figur 1) — diese 
-Grenze besser anzunahern erlauben. 
_ Der Widerstand der Basis-Emitterstrecke und der Sperrstrom J, des Kollek- 
tors sind stark temperaturabhangig. Ihre Werte sollten deshalb keinen nennens- 
werten Einfluss auf das Verhalten der Schaltung besitzen. Dies wird so erreicht, 
dass im Stromkreise des Emitters ein geniigend gross dimensionierter Seriewider- 
‘stand Re eingeschaltet wird und ferner der Strom durch den Spannungsteiler 
Rc i“ R, — R, gross gegen [,, bemessen wird. 
Figur 2 zeigt das Schema der Ausgangsstufe. Zur Registrierung dient ein 
Sodeco-Zahlwerk von 12 Q Widerstand. Der dazu notwendige Strom von etwa 


0c70 0C15 
4 ‘ 


alles 


+ 
Figur 2 
Endstufe mit Leistungstransistor und Zahlwerk. 


150 mA wird vom Leistungstransistor OC 15 geliefert. Um die letzte Untersetzer- 
stufe durch den Steuerstrom nicht tibermassig zu belasten, erfolgt die Kopplung 
iiber zwei Kollektorbasisstufen (Kathodefolgestufen). 

Die Betriebsdaten des Untersetzers sind der folgenden Zusammenstellung zu 


entnehmen: 


Transistoren OC 70 und 71 SB-100 
Pernt ew 22 Nh ibe xs! SEyt po 400 kHz 30 MHz 
Batteriespannung .-.-.-+--: : 655) V1 4,5 V 
Batteriestrom . . . » =» + + + 2,6mA 0,9mA 
Eingangsimpulse. .... +--+: ; 1S 11 Vv 
Ausgangsimpulse. . -- +++; ZOEN’ AeA 
Auflésungsvermogen 

2 Cg ath nm eo REZ 5 MHz 
BO Cote. eat eu AO? bite 5 MHz 
AU Noe fa culvaacey Calin 95 kHz 5) Marz 
LO, ee | Puce cy 80 kHz 5 MHz 


Zwischen statistischen und regelmassigen Eingangsimpulsen konnte kein Unter- 
schied im Auflésungsvermogen festgestellt werden. Anderungen der Batterie- 


spannung sind in grossen Grenzen zulassig. 
Wir danken der Stiftung Hasler- Werke, durch deren finanzielle Unterstutzung 


die vorliegende Arbeit zustande gekommen ist. 


as , ~ 
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Summary 


A simple circuit for a scale of two using junction transistors is described. 
A resolution up to 5 Mc/s has been achieved. 
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Elektromagnetische Wellen. Von Kart Witty WaGNER (Birkhauser — 
Verlag, Basel 1953). 257 S., 185 Fig., 22 Tab.; Lw. sFr. 33.30, br. sFr. 29.10. 

Nach langer, fruchtbarer und erfolgreicher wissenschaftlicher Tatigkeit, welche 
vor allem den theoretischen Grundlagen der elektrischen Nachrichtentechnik ge- 
widmet war, schloss KARL WILLY WAGNER im Jahre 1953, nur kurze Zeit nach © 
Vollendung seines 70. Altersjahres, die Augen fiir immer. Sein Buch Elektvo- 
magnetische Wellen ist eine Frucht der Arbeit seiner letzten Jahre. Es konnte 
noch vollkommen abgeschlossen und in mustergiiltiger Ausstattung in der Reihe | 
der Lehr- und Handbiicher der Ingenieurwissenschaften des Birkhauser Verlages 
in Basel erscheinen. Sein Zweck besteht vor allem in einer Einfiihrung in die | 
Theorie der elektromagnetischen Wellen als Grundlage fiir ihre Anwendungen | 
in der elektrischen Ubertragungstechnik. Es ist daher Wert auf eine sehr klare, | 
nicht zu knappe Darstellung gelegt, welche durch viele sorgfaltig ausgewdhlte 
Figuren und Tabellen unterstiitzt wird und die mathematischen Kenntnisse nicht 
tiberschreitet, die ein Ingenieur an einer Technischen Hochschule zu erwerben | 
pflegt. Den Bediirfnissen des Kabeltechnikers (Drahtwellen, Wellen in Leitungs- — 
biindeln) oder des Starkstromtechnikers (Wanderwellen) wird ebenso Rechnung 
getragen wie denjenigen des Radiotechnikers (Wellen in Hohlleitern, Dipole, 
elektrische Wellen in der Ionosphare). Die eingestreuten Anwendungen und Bei- : 
spiele, die oft bis zur zahlenmassigen Berechnung gehen, tragen zur Belebung des — 
Stoffes und zur Erleichterung des Verstandnisses erheblich bei. Man spiirt beim 
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Studium des Buches, dass es nicht nur von einem begabten Theoretiker verfasst 
wurde, sondern auch von einem durch lange, enge Verbindung mit der Praxis 
erfahrenen Fachmann. Die besten Empfehlungen des Referenten begleiten das 

Werk. F. Tank 


Neutron Optics. Von D. J. HUGHES (Interscience Publishers, Inc., New York 
1954). VII + 136 S., 34 Fig.; $2.50. 

Das Buch von D. J. HuGues befasst sich sowohl mit den physikalischen 
Prinzipien als auch mit den Anwendungen der Neutronenoptik. Dabei wird in 
erster Linie auf das Verhalten und die Benutzung von thermischen Neutronen 
eingegangen, die in Uranreaktoren in extrem hohen Dichten erzeugt werden 
kénnen. Das Buch ist wie folgt gegliedert: Das erste Kapitel umfasst die Grund- 
lagen der Neutronenoptik, wobei auf Verwandtes und Gegensatzliches zur ge- 
wohnlichen Optik hingewiesen wird. Eine kurze Ubersicht der experimentellen 
Methoden ist im folgenden Kapitel enthalten. Die restlichen Kapitel sind den 
Anwendungen gewidmet: Wechselwirkungen der Neutronen mit Kernen, Be- 
stimmung von Kristallstrukturen und Neutronenstreuung an magnetischen Ma- | 
terialien. Das letztgenannte Gebiet umfasst die pragnantesten Anwendungen der 
Neutronenoptik. In dieser Domane haben Neutronen gegentiber elektromagneti- 
schen Wellen den Vorzug, dass sie Wechselwirkungen mit den magnetischen 
Momenten der Atome zeigen. 

Das Buch von HuGuHEs vermittelt einen klaren und iibersichtlichen Einblick 
in ein in den letzten Jahren entstandenes Wissensgebiet. Durch Weglassen von 
experimentellen und theoretischen Einzelheiten wird ein vorziiglicher Uberblick 
gegeben. P. Huber 


Praktische Mathematik fiir Ingenieure und Physiker. Von R. ZURMUHL 
(Springer, Berlin 1953). XI + 481 S., 114 Abb.; DM 28.50. 

Dieses als Erginzung der mathematischen Grundvorlesung gedachte Buch 
ist eine Einfiihrung in die numerischen Methoden der Mathematik, unter beson- 
derer Beriicksichtigung der (elektrischen) Rechenmaschinen. Hauptgegenstande 
sind: Auflésung von Gleichungen, insbesondere von algebraischen und von Syste- 
men linearer, Matrizen-Eigenwertprobleme, Interpolation und Integration, Aus- 
gleichsrechnung, Reihenentwicklung von Funktionen und schliesslich die An- 
fangs-, Rand- und Eigenwertprobleme gewohnlicher Differentialgleichungen. In 
jedem Kapitel wird mit der mathematischen Theorie von den Grundlagen her 
begonnen, so dass das Buch schon in den ersten Semestern gelesen werden kann. 
Einige spezielle Fragen, wie zum Beispiel die Auflosung von Gleichungen dritten 
und vierten Grades, sind recht breit behandelt; dagegen wird mancher Leser, 
vor allem der Ingenieur in der Praxis, bedauern, dass wichtige Problemkreise, 
wie partielle Differentialgleichungen. oder Integralgleichungen, nicht mehr zur 
Behandlung kommen. 

Beispiele und Rechenschemata sind nach Art einer Vorlesung in den ausfithr- 
lichen Text eingeflochten. Das Buch eignet sich daher weniger als Nachschlage- 
werk als zum Beispiel dasjenige von MitneE (Numerical Calculus | Princeton 
University Press, Princeton 1949], 393 Seiten), das leider nicht zitiert wird, wie 
denn iiberhaupt das angelsachsische Schrifttum etwas stiefmiitterlich behan- 
delt ist. 

Im iibrigen aber ist das Buch sehr sorgfaltig geschrieben, enthalt viele wert- 
volle Ratschlage fiir das numerische Rechnen und wird als Lehrbuch vor allem 
fiir den. Anfanger gute Dienste leisten. H. J. Maehly 


S.; $4.00. 
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i: briss der Theorie (Zusammen- 
_ hange mit der hypergeometrischen Verteilung und der unvollstandigen Beta-— 
ee funktion), eine Anleitung zum Gebrauch der Tafeln mit Ablesebeispielen, eine | 
4 Formel zur Interpolation in m und # und ein kleines Literaturverzeichnis. ‘< 
- Das vorliegende Werk bildet die Fortsetzung einer analogen Tafel fiir 7 = _ 
_ 2 (1) 49: Tables of the Binomial Probability Distribution, Applied Mathematics 
Series 6 (U.S. Government Printing Office, Washington 1950). Beide Tafeln wer- 
den. bei der Auswertung statistischer Qualitatskontrollen gute Dienste leisten. 
nt ee ee : H. J. Maehly 
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Flow Properties of Disperse Systems. Herausgegeben von J. J. Hermans. 
_ (North-Holland Publishing Company, Amsterdam 1953), x1 + 445 S., 162 Figs; 3 
Seth 35.— (70's)s. ae . 
’ In ay von J. M. Burcers, J. J. HERMANS und G. W. Scott BLAIR verOffent- — 
lichten Biicherreihe ‘‘ Deformation and Flow — Monographs on the Rheological — 
_» Behaviour of Natural and Synthetic Products’? nimmt der von J. J. HERMANS © 
 herausgegebene Band Flow Properties of Disperse Systems eine markante Stellung — 
ein. Darin wird erstmals versucht, die gemeinsamen Merkmale im str6mungs- _ 
_. technischen Verhalten disperser Systeme trotz ihrer mannigfaltigen Erscheinungs- 
_  formen herauszuarbeiten und zu einer allgemeinen, zweckdienlichen Darstellung 
~ zu bringen. Der Erfolg ist in Anbetracht der Vielgestaltigkeit des behandelten — 
__ Stoffes erstaunlich, wobei dieser die Strémungscharakteristiken aller Formen 
_ disperser Systeme umfasst, von den Suspensionen, Emulsionen und hochpoly- — 
_° + meren Lésungen bis zu den Schaumen, Rauchen und Pulvern. Je nach Gegen- | 
_ stand ist die Behandlung des Stoffes verschiedenen F. achbearbeitern zugewiesen. — 
_- Bemerkenswert ist, dass trotz der Betonung der gemeinsamen Merkmale im str6- 
-. mungstechnischen Verhalten, die kennzeichnenden Eigenschaften der dispersen _ 
_. Partikel und ihr individueller Einfluss auf die Str6mungserscheinungen ebenfalls — 


.. zur gentigenden Darstellung gelangen. . 

Pb. Das Buch umfasst die folgenden Aufsditze: I. R. Roscor (Newcastle upon 
—.* Tyne), Suspensions; E. G. RICHARDSON (Newcastle upon Tyne), Emulsions: 
-* . J. J. Hermans (Leiden), Gels; CH. SADRON (Strasbourg), Dilute Solutions of 


ww 


_ Impenetrable Rigid Particles; J. J. HERMANS (Leiden), Dilute Solutions of Flexible 
» Chain Molecules; E. G. RicHaRDsON (Newcastle upon Tyne), Liquid Sprays; 
—. H.L. Green (Porton Down, Wilts.), Atomization of Liquids; RAPHAEL MaTALON. 
Rix’ _. (London), Foams; H. L. GREEN (Porton Down, Wilts.), Smoke; Mrs. B. S. Neu- 
MANN (Redhill, Surrey), Powders. : 

Den Schluss des Buches bildet ein ausfiihrliches Sach- und Autorenregister. — 

Das Werk lasst berechtigterweise erwarten, dass es schon in naher Zukunft zur 

Ausriistung aller Hochschul- und Industrielaboratorien der Physik, Chemie und — 

_., . Biologie gehéren wird, in denen das stromungstechnische Verhalten disperser 

-. Systeme studiert wird. R. Sanger 
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ELEKTRONENMIKROSKOPE 
MIKROTOME 
ELEKTRONENDIFFRAKTOGRAPHEN 
HOCHSPANNUNGSOSZILLOGRAPHEN 


Elementare Differentialgeometrie 


Von Dr. Woircanc Haack, Professor an der Technischen Hochschule Berlin. 242 Seiten 
mit 25 Abb. Preis gebunden Fr. 22.- (DM 22.—), brosch. Fr. 18.70 (DM 18.70). Sammlung 
«Lehrbiicher und Monographien aus dem Gebiete der exakten Wissenschaften» — Mathe- 
matische Reihe, Band 20. 

Inhalt: Vektoren und Bewegungsinvarianten — Kurven im Raum — Einparametrige 
Dreibeinscharen. Streifen — Elemente der Flachentheorie — Abbildungen zweier Flachen 
aufeinander — Ableitungsgleichungen und Integrierbarkeitsbedingungen — Geometrie 
auf der Flache — Minimalflachen — Zweidimensionale Dreibeinmannigfaltigkeiten — 
Existenzsatze der Flachentheorie. 


Weyl 
Symmetrie 


Von Dr. HERMANN WEYL, Professor an der Eidg. Technischen Hochschule in Ziirich und 
am Institute for Advanced Study in Princeton. Aus dem Englischen tibersetzt von Prof. 
Dr. Lutu BecuToisHeEm, University of Redlands/Cal. (USA.). 157 Seiten mit 72 Figuren. 
Preis gebunden Fr. 17.70 (DM 17.70). Sammlung «Wissenschaft und Kultur» — Band 11. 
Inhalt: 1. Bilaterale Symmetrie — 2. Translatorische, rotative und verwandte Symme- 
frien — 3. Ornamentale Symmetrie — 4. Die allgemeine mathematische Idee der Symme- 
trie — 5. Zwei Anhange tiber endliche Rotationsgruppen in drei Dimensionen. 

Wo immer der Mensch Ordnung, Schénheit und Vollkommenheit zu begreifen oder zu 
schaffen versucht hat, war Symmetrie ihm ein leitendes Prinzip. Dieses mit zahlreichen 
Abbildungen geschmiickte «Handbiichlein der Symmetrie», das in deutscher Ubersetzung 
vier an der Universitat Princeton 1951 gehaltene Vorlesungen (Vanuxem Lectures) 
wiedergibt, sucht auf kleinem Raum allen wesentlichen Seiten dieses faszinierenden 
Gegenstandes gerecht zu werden. 
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a Schlafli 
| ~ Gesammelte mathematische Abhandlungen 


In drei Banden. Herausgegeben im Auftrag der Schweizerischen Mathematischen Gesell- 
_ schaft vom Komitee Steiner-Schlafli der Schweizerischen Naturforschenden Gesellschaft, 

erganzt und erlautert durch Nachworte der Herausgeber. 

Band III: 402 Seiten. Preis gebunden Fr. 59.30 (DM 59.30). : 

_«Die Schweizer Kollegen sind zu begliickwiinschen, dass es ihnen gelungen ist, ihrem > 
hee grossen Landsmann mit dieser Gesamtausgabe seiner mathematischen Schriften ein wiir- 
yy diges Denkmal zu schaffen. Die Ausstattung und der Druck des Buches sind hervor- 
ragend» (H. Gortler, Archiv der Mathematik). 

“The elegant volume under review, profusely annotated by Burckhardt, Hadwiger, and 
Kollros, is a fitting tribute to the genius of this great Swiss mathematician” (H. S.M. - 
Coxeter, Mathematical Reviews). 

Mit diesem Band ist das dreibandige Werk, das erstmals samtliche Schriften von Bedeu- 
tung des grossen Mathematikers enthalt, abgeschlossen. 


Speiser 
Die geistige Arbeit 


Von Dr. ANDREAS SpEIsER, Professor an der Universitat Basel. 216 Seiten. Preis gebun- 
den Fr. 19.25 (DM 19.25), Sammlung «Wissenschaft und Kultur» — Band 9. 

Inhalt: Die geistige Arbeit — Die raumliche Deutung der Aussenwelt — Die platonische 
Lehre vom Unbekannten Gott und die christliche Trinitat — Platos Ideenlehre und die 
Mathematik — Wissenschaft und Glaube — Der Erlésungsbegriff bei Plotin — Die mathe- 
matische Betrachtung der Kunst — Die Gottesbeweise von Plato, Anselmus und Leibniz — 
Die philosophisch-naturwissenschaftliche Fakultat — Geist und Mathematik — Die Grund- 
lagen der Mathematik von Plato bis Fichte — Uber die Freiheit — Die mathematische 
Denkweise in den Geisteswissenschaften und in der Kunst. 
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ee Selecta Hermann Weyl 


Herausgegeben zum siebzigsten Geburtstag von Prof. Dr. HERMANN WEYL von der Eidg. 
Technischen Hochschule in Ziirich und vom Institute for Advanced Study in Princeton. 
Etwa 600 Seiten. Preis gebunden Fr. 48.90 (DM 48.90). 

Diese Sammlung enthalt 19 Abhandlungen von Hermann Wey] aus den Jahren 1910 bis 
1952; sie gibt einen Querschnitt durch das mathematische Schaffen des grossen Gelehrten. 
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